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Abstract 



In the heavy-quark hmit, the hadronic matrix elements entering nonleptonic B- 
meson decays into two light mesons can be calculated from first principles including 
"nonfactorizable" strong-interaction corrections. The B — > ttK, tttt decay ampli- 
tudes are computed including electroweak penguin contributions, SU(3) violation 
in the light-cone distribution amplitudes, and an estimate of power corrections 
from chirally-enhanced terms and annihilation graphs. The results are then used 
to reduce the theoretical uncertainties in determinations of the weak phases 7 and 
a. In that way, new constraints in the (p, ff) plane are derived. Predictions for the 
B ttK, tttt branching ratios and CP asymmetries are also presented. A good 
global fit to the (in part preliminary) experimental data on the branching fractions 
is obtained without taking recourse to phenomenological models. 



1 Introduction 



The study of nonleptonic two-body decays of B mesons is of primary importance for 
the exploration of CP violation and the determination of the flavour parameters of the 
Standard Model. Because of the interference of several competing amplitudes, these 
processes allow for the presence of different weak and strong-interaction phases, which 
play a crucial role for CP violation. In the Standard Model, all CP-violating observables 
are related to the complex phase of the quark mixing matrix, which in turn implies 
nontrivial angles in the "unitarity triangle" VudV*^, + VcdV*^, + VtdV^l = 0. With the 
standard choice of phase conventions, one defines the weak phases (3 = —SiTg{Vtd) and 
7 = arg(V^*^), as well as a = 180° — [3 — •y. In the Standard Model, sin 2/3 can be 
extracted in a theoretically clean way by measuring the time-dependent rates for the 
decays 5°, 5° J/ipKs- The measurement of 7 (or a) is more difficult, since it 
requires controlling the hadronic dynamics in nonleptonic B decays. 

A promising strategy for the determination of 7 is based on rate measurements for the 
charged modes B^ — >■ (irK)^ and B^ — > tt^tt^ Hadronic uncertainties in this method 
can be reduced to a minimum by exploiting the structure of the effective weak Hamilto- 
nian and using isospin and SU(3) flavour symmetries 0, |^, |, ^. If only measurements 
of CP-averaged branching ratios are available, it is still possible to derive bounds on 7 
0, 1^, or to determine 7 under the assumption of only a moderate strong-interaction 
phase between penguin and tree amplitudes. A different strategy for extracting 7 
uses SU(3)-symmetry relations between the various contributions to the time-dependent 
BdiBfi TT+TT^ and Bg^B^ K^K~ decay amplitudes 0. The main theoretical limi- 
tation of these methods is in the accuracy with which the effects of SU(3) breaking can 
be estimated. 

The angle a can be determined from the time- dependent CP asymmetry in the decays 
B^,B^ TT+TT^, if the subdominant penguin contribution to the decay amplitudes can 
be subtracted in some way. This can be achieved by an isospin analysis f^; however, 
in practice that route is extremely challenging due to the difficulty in measuring the 
very small B^,B^ ir^ir^ branching ratios. One must therefore rely on dynamical input 
for the penguin-to-tree ratio ||^, 0. Alternatively, a can be measured in related 
decays such a.s B p it, for which the penguin contribution can be eliminated using a 
time-dependent analysis of the B tt+tt^tt'^ Dalitz plot |r2| . 

Most of the above-mentioned determinations of 7 or a have theoretical limitations, 
which would be reduced if some degree of theoretical control over two-body nonleptonic 
B decays could be attained. However, in the past this has proven to be a very difficult 
problem. Even advanced methods such as lattice gauge theory, QCD sum rules, or the 
large-Nc expansion have little to say about the QCD dynamics relevant to hadronic 
B decays. In recent work, we have developed a systematic approach to this problem. 
It is based on the observation that, in the heavy-quark limit ^ Aqcd, a rigorous 
QCD factorization formula holds for the two-body decays B M1M2, if the "emission 
particle" M2 (the meson not obtaining the spectator quark from the B meson) is a 
light meson. (It has been argued that perhaps the large- iVc limit may be more relevant 
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to factorization than the heavy-quark hmit |T^. However, the dramatic decrease of 
"nonfactorizable" effects seen when comparing K ^ tttt, D ^ ttK and B — >■ ttK decays 
shows that the heavy-quark hmit is of crucial importance.) We have previously applied 
the factorization formula to i? ^ vrvr decays and obtained results for the decay amplitudes 
at next-to- leading order in and to leading power in AqcB/^b fl^- '^^^ conceptual 



foundations of our approach have been discussed in detail in [|T^, which focused on the 
simpler case of decays into heavy-light final states (where Mi is a charm meson). In 
the present work, the QCD factorization formula is applied to the general case of B 
decays into a pair of light, flavour-nonsinglet pseudoscalar mesons. (Preliminary results 



of this analysis have been presented in |jT6[.) The present analysis contains three new 



theoretical ingredients in addition to a much more detailed phenomenological analysis: 

1. Matrix elements of electroweak penguin operators are included, which play an im- 
portant role in charmless decays based on 6 — > sqq transitions. This is a straightforward 
extension of our previous analysis. However, a sensible implementation of QCD correc- 
tions to electroweak penguin matrix elements implies that one departs from the usual 
renormalizat ion-group counting, in which the initial conditions for the electroweak pen- 
guin operators at the scale fi = Mw are treated as a next-to-leading order effect. 

2. Hard-scattering kernels are derived for general, asymmetric meson light-cone dis- 
tribution amplitudes. This is important for addressing the question of nonfactorizable 
SU(3)-breaking corrections, since the distribution amplitudes of strange mesons are, in 
general, not symmetric with respect to the quark and antiquark momenta. 

3. The leading power corrections to the heavy-quark limit are estimated by analyzing 
"chirally-enhanced" power corrections related to certain twist-3 distribution amplitudes 
for pseudoscalar mesons. We also discuss potentially large power corrections arising from 



annihilation topologies, noted first in This is essential for controlling the theoretical 
uncertainties of our approach. 

The second and third items have not been considered in previous generalizations of the 
results of to the case of -B — nK decays [1T8| , [T9| , |20| , Chirally-enhanced power 



corrections were discussed in [|^, |23[ , but we disagree with some of the results obtained 
by these authors. 

The QCD factorization approach provides us with model-independent predictions 
for the decays amplitudes including, in particular, their strong-interaction phases. The 
same method can also be applied to other charmless decays, such as vector-pseudoscalar 
||T9| or vector-vector modes. Our main focus here is on the development of the new 
conceptual aspects of the approach that are important for a comprehensive phenomeno- 
logical analysis. This includes a detailed discussion of various sources of potentially large 
power corrections to the heavy-quark limit. We then perform a comprehensive study of 
CP-averaged branching fractions and direct CP asymmetries in decays to t^K and tttt 
final states, including a detailed discussion of the theoretical uncertainties from all in- 
puts to the QCD factorization approach. In many of the phenomenological applications 
discussed in this work the dynamical information obtained using the QCD factorization 
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formalism is used in a "minimal way", to reduce the hadronic uncertainties in methods 
that are theoretically clean up to "nonfactorizable" SU(3)-breaking effects. Most im- 
portantly, these strategies do not suffer from uncertainties related to weak annihilation 
contributions. In this way, it is possible to reduce the hadronic uncertainties in the 
strategies for determining 7 from nK, Tin decays proposed in [|I], Q to the level of 

"nonfactorizable" corrections that simultaneously violate SU(3) symmetry and are power 
suppressed in the heavy-quark limit. These corrections are parametrically suppressed by 
the product of three small quantities: l/Nc, {nis — rrid) / Aqqb, and Aqcd/"^;)- As a conse- 
quence, we argue it will eventually be possible to determine 7 with a theoretical accuracy 
of about 10° (unless 7 is much different from its expected value in the Standard Model). 
More accurately, the strategies discussed here can constrain the Wolfenstein parameters 
p and f] with accuracies similar those obtained from the standard global fit to \Vub\, ^k, 
and B-B mixing. Given the theoretical input discussed in this paper even the present, 
preliminary data on the rare hadronic decays exclude at 95% confidence level half of the 
parameter space obtained from the standard fit. 

The QCD factorization approach provides a complete theoretical description of all 
B —>■ PP decay amplitudes. This allows for a large variety of predictions, which go far 
beyond those explored in the present work. In the future, this will offer the possibility 
of several nontrivial experimental tests of the factorization formula. A more detailed 
discussion of these predictions will be presented elsewhere. 

The remainder of this paper is organized as follows: In Section 0, we collect some 
basic formulae and express the B ttK, tttt decay amplitudes in terms of parameters 
a, and bi appearing in the effective, factorized transition operators for these decays 
(including weak annihilation contributions). Section ^ contains the technical details of 
the calculations based on the factorization formula, a discussion of annihilation effects, 
and a compilation of the relevant formulae for the numerical evaluation of our results. 
Readers not interested in the technical aspects of our work can proceed directly to 
Sections | and ^, where we present numerical values for the amplitude parameters a, 
and bi (Section ^ and discuss phenomenological applications of our results (Section |^). 
Specifically, we consider strategies to bound and determine the weak phase 7 and to 
extract sin 2a from mixing- induced CP violation in B tt+tt" decay. We also present 
predictions for CP-averaged branching fractions and CP asymmetries, and perform a 
global fit in the (p, fj) plane to all measured B ttK, tttt branching fractions. A critical 
comparison of our formalism with other theoretical approaches to hadronic B decays is 
performed in Section |. 

2 Parameterizations of the decay amplitudes 

The effective weak Hamiltonian for charmless hadronic B decays consists of a sum of 
local operators Qi multiplied by short- distance coefficients Cj and products of elements 
of the quark mixing matrix, Ap = VphV*^ or = VphV*^. Below we will focus on B — > nK 
decays to be specific; however, with obvious substitutions a similar discussion holds for 
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all other B decays into two light, flavour- nonsinglet pseudoscalar mesons. Using the 
unitarity relation — Aj = + Ac, we write 

^cff = % E QI + C2QI+ E + C77 Ql, + Cs, Qsg) + h.C. , (1) 

V ^ p=u,c ^ j=3,...,10 ^ 

where Q\ 2 are the left-handed current-current operators arising from VF-boson exchange, 
Qz,...fi and (57,...,io are QCD and electroweak penguin operators, and Qy^ and Q^g are the 
electromagnetic and chromomagnetic dipole operators. They are given by 

Ql = {ph)v-A{sp)v~A , Ql = ipibj)v~AiSjPi)v-A , 

Qs = isb)v-A E? i^Q)v-A , Q4 = iSibj)v-A E? iW-^v-A , 

Q5 = (s6)v-A E? {^(l)v+A , Qe = E? , 

Q7 = {sh)v-A E? ie<7(gg)y+A , Qs = {sihj)y_A Ee ie<?(gigi)y+A , 

Qg = (s6)v_A E? ieg(f7?)y-A , QlO = (S»&j)y-A E? l^qi.mi)v-A , 

= ^m6sa^,(l + 75)i"'^"&, Qsg = ^ m,m^,{l + ^^)G^''h , (2) 

where (^ig2)y±A = ^i7aj(1 ±75)'?2, J are colour indices, Cq are the electric charges of the 
quarks in units of |e|, and a summation over q = u,d, s, c, b is implied. (The definition of 
the dipole operators Qy-y and Qsg corresponds to the sign convention iD^^ = id^ + QgA^ta 
for the gauge- covariant derivative.) The Wilson coefficients are calculated at a high scale 
/i ~ My/ and evolved down to a characteristic scale /i ~ using next-to-leading order 
renormalization-group equations. The essential problem obstructing the calculation of 
nonleptonic decay amplitudes resides in the evaluation of the hadronic matrix elements 
of the local operators contained in the effective Hamiltonian. 

Applying the QCD factorization formula and neglecting power-suppressed effects, the 



matrix elements of the effective weak Hamiltonian can be written in the form 14, 15 



where 



(vririKffl^) = %T.K + T;'^^\B) , (3) 

V ^ p=u,c 



Tp = aii^nK) Spu {ub)v-A ® isu)v-A 
+ a2{'nK) 6pu {sb)v-A ® {uu)v-a 
+ a^{TxK) Eg {sb)v-A ® {qq)v-A 
al{7iK) E<? {qb)v-A ® {sq)v-A 
+ a^i^TiK) Eg {sb)v-A ® {qq)v+A 
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+ al{7cK) {-2){qb)s^P ® {sq)s+P 
+ a-ji^-nK) {sh)v-A ® leq{qq)v+A 
+ al{TiK) Y,, {-2){qh)s-P ® |e,(sg)5+P 
+ agiT^K) {sh)v-A ® leq{qq)v-A 

+ al^i^nK) Y,q {qb)v-A ® f eg(sg)y_A • (4) 

Here {qiq2)s±p = i 75)q'2, and a summation over q = u,d is implied. The symbol 
® indicates that the matrix elements of the operators in Tp are to be evaluated in the 
factorized form {7iK\ji ® j2\B) = {7i\ji\B) {K\j2\0) or {K\ji\B) (7r|j2|0), as appropriate. 
"Nonfactorizable" corrections are, by definition, included in the coefficients Oj. The 
matrix elements for B mesons (i.e., mesons containing a 6-antiquark) are obtained from 
(|) by CP conjugation. A corresponding result, with obvious substitutions, holds for 
other decays such as -B ^ tttt. 

The term 7^^"^° in arises from weak annihilation contributions and introduces a 



set of coefficients bi{7rK), which we shall define and discuss in detail in Section ^ 
Annihilation contributions are suppressed by a power of Aqcd/"'^?) and not calculable 
within the QCD factorization approach. Nevertheless, we will include the coefficients bi 
in the amplitude parameterizations in this section. 

The coefficients Oj multiplying products of vector or axial-vector currents are renor- 
malization-scheme invariant, as are the hadronic matrix elements of these currents. For 
the coefficients and as a scheme dependence remains, which exactly compensates 
the scheme dependence of the hadronic matrix elements of the scalar or pseudoscalar 
densities associated with these coefficients. These matrix elements are power suppressed 
by the ratio 

^i^(^) = '^Hhi ^ (5) 

^ nibifi) {mq{n) + msifi)) ' 

which is formally of order Aqqd /rrih but numerically close to unity. In the following we 
shall use the same notation for charged {q = u) and neutral kaons (g = ci), since the 
difference is tiny. A corresponding ratio 

2m? 

^^^^^ mbifi) (niuifi) +md{fi)) ^^"^ 

appears in the discussion of the B ^ titt decay amplitudes. For a phenomenological 
analysis of nonleptonic B decays it is necessary to estimate these "chirally-enhanced" 
corrections despite the fact that they are formally power suppressed. A detailed discus- 
sion of these corrections will be presented in Section |3.2| , and their numerical importance 
will be investigated in Section |^. 

In terms of the parameters a,, the B —>■ nK decay amplitudes (without annihilation 
contributions) are expressed as 



K 
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-V2AiB- - 



T^-\ 



■kK 



+ 



Am O2 + -^p 2 + 0.9) 



Aktt } 



■K+K 
.0 



A„ ai + Ap {al + afg) + Ap (a^ + a; 



A{B' 



(7) 



Here Ap = VpbV*g, = aiijiK), and a summation over p = m, c is implicitly understood 
in expressions like Apof. The last relation follows from isospin symmetry. The CP- 
conjugate decay amplitudes are obtained from the above by replacing Ap A*. We have 
defined the factorized matrix elements 
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v2 



^k) Ik , 



•-KJ ^ 



(8) 



where FQ^^{q'^) are semileptonic form factors. Weak annihilation effects contribute 
further terms to the decay amplitudes, which can be parameterized as 



where 
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A„62 + (A„ + Ae)(63 + &3 ) 
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) •^a.nniB > 

K = fsfirfK ■ 



2 



B 



nK 



(9) 
(10) 



The coefficients bi = bi{TTK) will be defined in Section p.5| , but we may note here that 
bi^2 are related to the current-current operators Qi and Q2 in the effective Hamiltonian 
(H), and 63,4 (b^Y) related to QCD (electroweak) penguin operators. 
The i? — s> TTTT decay amplitudes are given by 



TT TT 



At, 



A^ fli + Ap {al + a?o) + A^ (ag + a 
3 

A^(ai + 02) + - Ap(-a7 + + ag + a- 



p \ 
10) 



At, 



TT 



(11) 
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where now A' = VpbV*^, ai = aii^mr), and 



A.. = (m| - ml) F,^--{ml) U . (12) 



The additional annihilation contributions are 
-72 Ann(5- ^ 7r-7r°) =0, 

Ann(5° ^ VrOvrO) = -Ann(5° ^ TT+TT-) , (13) 



where 6,- = 6, (tttt), and 



B^^ = ^^fBfl (14) 



Neglecting tiny mass corrections of order (m^r./f/wi^) 
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i?^^ = ^ ^ -^0 ^ ^ ^ . (15) 

These ratios will play an important role in the discussion of SU(3) violations in Section ^ 
The expressions collected above provide a complete description of the decay am- 
plitudes in terms of the parameters a, and bi for the various processes. They are the 
basis for most of the phenomenological applications discussed in this work. However, 
in the literature several alternative parameterizations of the B — > tiK decay amplitudes 
have been introduced, which are sometimes useful when considering CP asymmetries or 
ratios of branching fractions. We briefly elaborate on one such parameterization here, 
adopting the notations of 0]. The dominant contributions to the B —>■ nK decay am- 
plitudes come from QCD penguin operators. Because the corresponding operators in 
the effective weak Hamiltonian preserve isospin, this contribution is the same (up to 
trivial Clebsch-Gordon coefficients) for all decay modes. Isospin-violating contributions 
to the decay amplitudes are subdominant and arise from the current-current operators 
Qi and Q2 (so-called "tree" contributions), and from electroweak penguins. The latter 
are suppressed by a power of a/ whereas the former are suppressed by the ratio 

eKM e-'^ = hL = tan^^c Rb e"'^ , (16) 

Ac 

where 9c is the Cabibbo angle, 

R, = cot ^ = + (17) 



cb\ 



7 



is one of the sides of the unitarity triangle, and p and rj are the Wolfenstein parameters. 
A general parameterization of the decay amplitudes is 



A{B- 
-V2A{B- 



> -R-K 



0\ 



P 
P 



-«7 



£3/2 e 



qe" 



St e'^^ie-'^ 



qce 



(18) 



The amplitude A{B^ —>■ n^K^) is then de- 



(The phase e**^" was denoted — e*'' in [| 
termined by the isospin relation shown in the last line of (|^. The dominant penguin 
amplitude P is defined as the sum of all contributions to the B~ —>■ tt~K^ amplitude that 
are not proportional to e~^"' . This quantity cancels whenever one takes ratios of decay 
amplitudes, such as CP asymmetries or ratios of branching fractions. The parameters 
£3/2 and St measure the relative strength of tree and QCD penguin contributions, q and 
qc measure the relative strength of electroweak penguin and tree contributions, and Sa 
parameterizes a rescattering contribution to the B~ — >■ tt'K^ amplitude arising from 
up-quark penguin topologies. Moreover, 0, 0^, a;, o^c, and (pa are strong rescattering 
phases. (The strong- interaction phase of P is not observable and can be set to zero.) 
All parameters except q e*"^ receive weak annihilation contributions. 

It is straightforward to express the various amplitude parameters in terms of the 
parameters = ai{7rK) and bi = bi{iTK) defined earlier. We obtain 



P=Xr 



1 

r 



10 



+ r 



K 



1 . 

2 ^ 



K 



hf^) B 



■>ttK 



(19) 



for the leading penguin amplitude, and for the remaining parameters 

(ai + R^Ka'2) + iKo + r^al + Rnxjag - Q7)] 



£3/2 e 



icf) 



10 



3 uEW\ 



eaC 



-CRM 



ai + |«o + <<)-r^(62 + i&3 



« + r^O - i(a^o + r^a^) + r^(62 + &3 + 



+ rfa§)+rA(63 + 6fW) 



qe" 



no 



qce 



2eKM 
3 



(ai + R-nKa2) + |[aio + r^al, + R^xiag - a?)] 
+ rf ao - TAbl 



(20) 



2eKM a, + liafo + r^a^) - rA{b2 + | &f^) ' 
where Rt^k = Axn/AT^K is the ratio of the two factorized amplitudes given in (p!5D , and 
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nK 



ml Ft 



-(0) 



(21) 
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Our notation for amplitude ratios is such that the ratio Rt^k (denoted by a capital R) de- 
viates from 1 only by SU(3)-breaking corrections, whereas the ratios and (denoted 
by a lower-case r) are formally of order Aqco/'^fe in the heavy-quark limit. However, 
whereas ta ~ 0.003 is indeed very small, ^ 0.7 (at a scale ~ 1.45 GeV) is numer- 
ically large for real B mesons. Finally, note that the electroweak penguin coefficients 
O7,...,io could be safely neglected in all quantities other than ge*"^ and qce^'^'^, because 
they are tiny compared with the other coefficients ai_..._6- In (|20|), they are included only 
for completeness. (The systematics of including electroweak penguin contributions will 
be discussed in more detail later.) 

An important quantity affecting the determination of sin 2a from the time-dependent 
CP asymmetry in the decays -B°, 7r~^n~ is the ratio of penguin to tree amplitudes. 

In this case the tree contribution is no longer CKM suppressed, since 

^ = -R,e-'^ (22) 

is of order unity. The n'^ii^ decay amplitude in (p!T]) can be written as 

- ^(5° ^ TT+TT-) oc e-*^ + ^ , (23) 

where 

P.n 1 {at + r;a^) + (a^p + r^a^) + VAjh + 2h - |(&r - ^D] 

T^. Rb (ai + al + rjO + (a^^o + rja^) + r^[6i + h + 2h - |(6f^ - &f^)] ' ^ ^ 

and ta = Bt,t,/At,t,. In this case the electroweak penguin terms are very small, because 
they are not CKM enhanced with respect to the tree contribution. 

This concludes the discussion of parameterizations of the decay amplitudes. The 
following section is devoted to a detailed description of the QCD factorization formalism, 
the calculation of the parameters a, for the various nonleptonic decay amplitudes, and 
an estimation of the annihilation parameters hi. The reader mainly interested in the 
phenomenological applications of our results can proceed directly to Section ^ where 
we present numerical results for these parameters, which will be used in the subsequent 
analysis in Section ^ 



3 QCD factorization in S — > irK decays 

Based on the underlying physical principle of colour transparency (see |2B| for 

early discussions in the context of decays to heavy-light final states), supported by a 
detailed diagrammatic analysis of infrared cancellations at leading power in the heavy- 
quark expansion, we have shown in previous work that the complexity of the hadronic 
matrix elements governing energetic, two-body hadronic decays of B mesons simplifies 



greatly in the heavy-quark limit nib ^ Aqcd ||14|, [T5[- Consider B —>■ irK decays as an 
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Figure 1: Graphical representation of the factorization formula. Only one of 
the two form-factor terms in (^5]) is shown for simplicity. 



example. To leading power in Aqcd/'^6, but to all orders in perturbation theory, the 
matrix elements of the local operators Qi in the effective weak Hamiltonian in ([l|) obey 
the factorization formula 

{tiK\Q,\B) = Fff^^T^,, * fK^K + Fi^^Tl^^ * 

+ Tj'*fB^B*fK^K*U^.. (25) 

where are leading-twist light-cone distribution amplitudes, and the *-products imply 
an integration over the light-cone momentum fractions of the constituent quarks inside 
the mesons. A graphical representation of this result is shown in Figure Because 
the energetic, collinear light-quark pair that ultimately evolves into the emission particle 
at the "upper vertex" is created by a point-like source, soft gluon exchange between 
this pair and the other quarks in the decay is power suppressed in the heavy-quark 
limit (colour transparency). In other words, whereas the hadronic physics governing 
the semileptonic B — > Mi transition and the formation of the emission particle M2 is 
genuinely nonperturbative, "nonfactorizable" interactions connecting the two systems 
are dominated by hard gluon exchange. 

The hard-scattering kernels T/'^^ in (pST) are calculable in perturbation theory. T\ 



■ M,i 

starts at tree level and, at higher order in a^, contains "nonfactorizable" corrections from 
hard gluon exchange or light-quark loops (penguin topologies). Hard, "nonfactorizable" 
interactions involving the spectator quark are part of Tj^. The relevant Feynman dia- 
grams contributing to these kernels at next-to-leading are shown in Figure ^ Although 
individually these graphs contain infrared-sensitive regions at leading power, all soft and 
collinear divergences cancel in their sum, thus yielding a calculable short-distance con- 
tribution. Annihilation topologies are not included in (^) and Figure H, because they 
do not contribute at leading order in AQCD/''7ib. These power-suppressed contributions 



will be discussed separately in Section 3.5. 



We stress that the factorization formula does not imply that hadronic B decays are 
perturbative in nature. Dominant soft contributions to the decay amplitudes exist, which 
cannot be controlled in a hard-scattering approach. However, at leading power all these 
nonperturbative effects are contained in the semileptonic form factors and light-cone 
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A/ ^ \A 




o o 
o o 



Figure 2: Order ag corrections to the hard-scattering kernels T^j^ (first two 
rows) and T/^ (last row). In the case of Tlf-, the spectator quark does not 
participate in the hard interaction and is not drawn. The two lines directed 
upwards represent the quarks that make up one of the light mesons (the emis- 
sion particle) in the final state. 



distribution amplitudes. Once these quantities are given, the nonleptonic decay ampli- 
tudes can be derived using perturbative approximations to the hard-scattering kernels. 
This allows us to compute perturbative corrections to "naive factorization" estimates 
of nonleptonic amplitudes, which is crucial for obtaining results that are independent 
of the renormalization scheme adopted in the calculation of the effective weak Hamil- 
tonian. The hard-scattering kernels also contain imaginary parts, which determine the 
strong rescattering phases of the decay amplitudes. At leading power in Aqcd/"^6 these 
imaginary parts are of perturbative origin. 

In the remainder of this section, we discuss in detail the various issues to be addressed 



in the evaluation of the factorization formula. In Section a modified renormalization- 
group treatment of electroweak penguin effects is introduced, which is more appropriate 
than the standard scheme as far as applications to rare hadronic B decays are concerned. 
Each of the diagrams in Figure ^ contains a leading-power contribution relevant to (^) 
and power-suppressed terms, which do not factorize in general. An important class 
of such power-suppressed effects is related to certain higher-twist meson distribution 
amplitudes. These amplitudes are defined in Section |3^, and their leading, chirally- 



enhanced contributions to the nonleptonic decay amplitudes are evaluated. Section 
contains a compendium of the relevant formulae for the calculation of the parameters Oj. 
Section ^]5| is devoted to annihilation topologies and the definition of the parameters 6j. 



3.1 Wilson coefficients of electroweak penguin operators 

In the conventional treatment of the effective weak Hamiltonian the initial conditions 
for the electroweak penguin coefficients at the scale /i = M^r are considered a next-to- 
leading order effect, because they are proportional to the electroweak gauge coupling 
a (see [^, ^] for a detailed discussion). For our purposes, however, it is preferable 



11 



to deviate from this standard power counting and introduce a modified approximation 
scheme. The reason is that the electroweak penguin contributions in — >■ ttK decays 
and other rare processes based on 6 — > sqq transitions are important only because they 
compete with strongly CKM-suppressed tree topologies. Electroweak penguins and tree 
topologies together are responsible for the isospin-violating contributions to the decay 
amplitudes p9|]. Therefore, their effects are important even though they are suppressed 



with respect to the leading QCD penguin amplitude, which conserves isospin. The 
ratio of electroweak penguin to tree amplitudes scales like a/\^ ~ 1, where A = 0.22 
is the Wolfenstein parameter. Moreover, the dominant electroweak penguin effects are 
enhanced by a factor of {mt/MwY and 1/ ?,\v?6w- Hence it is not appropriate to count 
a as a small parameter in the renormalizat ion-group evolution, if the effect of interest is 
related to isospin breaking. 

We now describe a systematic modification of the usual leading and next-to-leading 
approximations, in which the dominant part of the electroweak penguin coefficients at 
the scale /i = Mt^ is treated as a leading-order effect. It is then consistent to include 
the QCD radiative corrections to the enhanced terms in the initial conditions for the 
electroweak penguin coefficients and, at the same time, the corrections of order a<j to 
the matrix elements of the electroweak penguin operators, which represent the next-to- 
leading order corrections to the hard-scattering kernels in the factorization formula. 

Using a compact matrix notation, the solution to the renormalization-group equation 
for the Wilson coefficients Ci, . . . , Cio in (||) can be written as 



C(/<) 



4:71 4:71 




C{Mw) . (26) 



The matrices Uq, J, Ro, and Ri depend on the ratio as{fi) / as{Mw) and on the anoma- 
lous dimensions and /^-function. At leading order, the evolution matrix reduces to 
Uq + {a/as)Ro. The remaining terms shown above are the next-to-leading corrections. 
We now expand the coefficients C{Mw) at the weak scale as 

CiM^) = + Ci'^ + £ ) + ^(^) + , (27) 

where superscripts indicate the order in the strong coupling constant as{Mw). The 
term proportional to a represents the electroweak contribution originating from photon- 
penguin, Z-penguin and box diagrams. We split this term into a contribution Ce con- 
taining all terms enhanced by the large top-quark mass and/or a factor of 1/ sin^6'iy, and 
a remainder Rg. As explained above, we treat Ce as a leading effect and hence include 
the first two terms in its expansion in powers of as{Mw)- The remainder Re is considered 
a next-to- leading effect, and so we only keep the first term in its perturbative expansion. 
Explicitly, the nonvanishing contributions to the initial conditions in the electroweak 
sector (z = 7, 8, 9, 10) are 

cSi = ^ , eg = ^ + [wBoixt) - 4Co{xt)] , (28) 

o o o sm '- -' 
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R^S = <9 = o Coixt) + - Doixt) - 4 , (29) 



and 

where Xt = m'f/M^ with rrit = mt{mt). The Inami-Lim functions Bq{x), Co{x) and 
Dq{x) can be found, e.g., in |]30[. Numerically, the remainder R'^'^ is indeed much 



smaller than C^^\ justifying our approximation scheme. Note that C^^^ is gauge and 
renormalization-scheme independent. The remainder R^^ is gauge-independent, but it 
carries the usual next-to-leading order scheme dependence of the electroweak coefficients. 
Explicit expressions for the QCD corrections contributing to (7^^^ have been obtained in 
||3T| . Using these results, we obtain the approximate expressions (valid for a high-energy 
matching scale fiw = Mw) 

eg ~ -29.60 28.52 

C«^0.94xr\ 

eg ~ -571.62x°-^«° + 566.40 

ego =^ -5.51 (30) 

In the conventional treatment C^^^ would be absent, while the sum (C^^^ + R^^^) would be 
the usual initial condition for the electroweak coefficients at the scale /i = M^, counted 
next-to- leading order effect. 

In addition to the modified counting scheme for powers of coupling constants, we make 
a further, numerically excellent approximation, which greatly simplifies the systematic 
evaluation of the Wilson coefficients. In essence, it amounts to neglecting QED effects in 
the calculation of the Wilson coefficients Ci, . . . , Cq. The values of these coefficients are 
obtained using the standard next-to-leading order approximation including only strong- 
interaction effects. At the same time, we neglect QED corrections to the matrix elements 
of the operators Qi, . . . ,Qq. (In fact, the virtual corrections of order a are infrared 
divergent and require the inclusion of photon bremsstrahlung contributions in order to 
obtain physical results. Our approximation scheme avoids this complication.) This 
treatment can be justified by noting that QED and electroweak contributions to the 
decay amplitudes in (|T8|) are only important if they contribute to the parameters ge*'^ 
and qce'^'^^'. From (0), it follows that the terms of order a contained in the coefficients 
0'7,...,io are enhanced by the prefactor l/exM- It is thus sufficient for all practical purposes 
to only include the order a corrections from the coefficients a7....,io- On the contrary, 
QED corrections to the other amplitude parameters can be safely neglected (i.e., the 
coefficients ai,...,6 do not contain terms proportional to a). Precisely this is achieved by 
our approximation scheme. 

At the technical level, the approximation described above can be explained in terms 
of the 10 X 10 anomalous- dimension matrix for the operators Qi, . . . ,Qio, written in 
block form as 



7 



^6x6 


-^6x4 







(31) 
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Table 1: Wilson coefficients Ci in the NDR scheme and based on our modified 
approximation scheme (see text). Input parameters are = 0.225 GeV, 
mt{mt) = 167 GeV, mfe(mfe) = 4.2 GeV, Mw = 80.4 GeV, a = 1/129, and 
sm'^Ow = 0-23. 



NLO 


Ci 


C2 




C4 


C5 


Ce 


fi = mb/2 


1.137 


-0.295 


0.021 


-0.051 


0.010 


-0.065 


fi = mb 


1.081 


-0.190 


0.014 


-0.036 


0.009 


-0.042 


fi = 2mb 


1.045 


-0.113 


0.009 


-0.025 


0.007 


-0.027 








Cg/a 


Cio/a 






fi = 'mb/2 


-0.024 


0.096 


-1.325 


0.331 







fi = rrib 


-0.011 


0.060 


-1.254 


0.223 








fi = 2mb 


0.011 


0.039 


-1.195 


0.144 








LO 


Ci 




C3 


C4 






fi = mb/2 


1.185 


-0.387 


0.018 


-0.038 


0.010 


-0.053 


fi = nib 


1.117 


-0.268 


0.012 


-0.027 


0.008 


-0.034 


fi = 2mb 


1.074 


-0.181 


0.008 


-0.019 


0.006 


-0.022 








Cg/a 


Cio/a 






fi = mb/2 


-0.012 


0.045 


-1.358 


0.418 


-0.364 


-0.169 


fi = mb 


-0.001 


0.029 


-1.276 


0.288 


-0.318 


-0.151 


fi = 2mb 


0.018 


0.019 


-1.212 


0.193 


-0.281 


-0.136 



We set C = 0, thereby neglecting the mixing of the electroweak penguin operators into 
the operators Qi, . . . , Qq, and ignore contributions of order a to A. At the same time, 
we drop the terms of order a in the matching conditions (p7|) for Ci, . . . ,Cq. We also 
omit terms of order a in D, which would yield second-order corrections in a. For the 
matrix B, we use the complete next-to- leading order result including terms of order a. 

Numerical results for the Wilson coefficients obtained at leading and next-to-leading 
order in our modified approximation scheme are given in Table Throughout this work, 
we use the "naive dimensional regularization" (NDR) scheme with anticommuting 75, 
as defined in |2^. The matrix elements of the dipole operators Qj^ and Qgg enter the 
decay amplitudes only at next-to-leading order. Consequently, the standard leading- 
logarithmic approximation is sufficient for the coefficients Cj^ and Cgg- In practice, it is 
advantageous to work with so-called "effective" coefficients, which in the NDR scheme 
are defined as C^^ = Cj^ — | C5 — Ce and C|J = Csg + C5. In the numerical analysis of 
(p6|) we consistently drop all terms of higher than next-to-leading order according to our 
modified counting scheme. Throughout we use the two-loop expression for the running 
coupling as{fi) evaluated with n/ = 5 light quark fiavours. 
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3.2 Meson distribution amplitudes and twist-3 projections 

Referring to the factorization formula shown graphically in Figure [l|, we denote by x 
the longitudinal momentum fraction of the constituent quark in the emission meson M2 
(the meson at the "upper vertex"), and by y the momentum fraction of the quark in 
the meson Mi. For a B meson decaying into two light mesons, we define light-cone 
distribution amphtudes by choosing the + direction along the decay path of the light 
emission particle and denote by ^ the light-cone momentum fraction of the light spectator 
antiquark. A massive pseudoscalar meson has two leading-twist light-cone distribution 
amplitudes [|15], but only one of them enters our results. This amplitude is called 



coincides with the function ^bi{0 defined in Section 2.3.3 of |T^. The meson 
distribution amplitudes are normalized to 1 once the decay constants are factored out as 
in For a light meson, we define the leading-twist amplitude in the usual way 
and assume that = 0(1) if both x and [1 — x) are of order unity, and $(x) = 0{x) 
for x — > (and similarly for x — > 1). For the B meson, almost all momentum is carried 
by the heavy quark, and hence $b(0 = 0{mb/ Aqcb) and ^ = 0{AQCT)/^b)- 

Higher-twist light-cone distribution amplitudes for the light mesons give power-sup- 
pressed contributions in the heavy-quark limit. However, as has been explained in Sec- 
tion ^ these can sometimes be large if they appear in conjunction with the chiral en- 
hancement factors r^(/i) defined in (^ and (|^). The corresponding terms are associated 
with twist-3 quark-antiquark distribution amplitudes and can be identified completely. 
Since the calculation of these contributions in momentum space is less straightforward 
than for the leading-twist contributions, we summarize the relevant projection operators 
below, following the discussion in [p^ . 

The relevant definitions of the light-cone distribution amplitudes of a light pseu- 
doscalar meson P in terms of bifocal operator matrix elements are ||3^ 



{P{pMz2h^l,q{z,)\0) = -ifpp^ ['dxe'^^^P-^'+'P-'^^^ 

•J 

JO 



X) 



{P{pMz2)a,,l,q{zi)\0) =tfpfip{p,z,-p,z,) j ' dx e'^'^-'^+'P-^^ , (32) 







where fp is the decay constant, and we have defined z = Z2 — zi and x = 1 — x. The 
parameter /ip = m?p/{mi + ^2), where mi^2 are the current quark masses of the meson 
constituents, is proportional to the chiral quark condensate. (This definition does not 
hold for the 7r° meson, in which case /x^o = m\ / {niu + md) as for the charged pions 
||T5| .) $(x) is the leading-twist (twist-2) distribution amplitude, whereas ^p{x) and 
have subleading twist (twist-3). All three distribution amplitudes are normalized 
to 1, as follows by taking the limit zi ^ Z2- The above definitions can be combined into 
the matrix 

{P{p)Mz2)qMm 
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We implicitly assume that the bilocal matrix elements are supplied with the appropriate 
path-ordered exponentials of gluon fields so as to make the definitions of the light-cone 
distribution amplitudes gauge invariant. (These exponentials are absent in light-cone 
gauge.) The distribution amplitudes depend on the renormalization scale /i. This scale 
dependence is compensated by higher-order corrections to the hard-scattering kernels, 
which however are beyond the accuracy of the present calculation. We thus suppress the 
argument /i in the distribution amplitudes, because it is irrelevant to our discussion. 

To obtain the corresponding projector of the quark-ant iquark amplitude in momen- 
tum space, the transverse components of the coordinate z must be taken into account. 
The collinear approximation can be taken only after the projection has been applied. 
We therefore assign momenta 

P P 

= xp^ + k^^ + ^ p", k^ = xp" - fc'f + _ ^ _ (34) 

2xp-p 2xp-p 

to the quark and antiquark in the light meson, where p is a light-like vector whose 3- 



components point into the opposite direction of p. Then the exponential in ( p3D becomes 
gi(fci-22+A:2-2i)^ (Meson mass effects are neglected, so that p and p can be considered as 
light-like.) The transverse components kj_ are defined with respect to the vectors p and 
p. Note that kf = k"^ = 0. In general, the projector (^) is part of a diagram expressed 
in configuration space. Transformation to momentum space is achieved by performing 
the integrations over Zi, which reduce to momentum-space 5-functions after substituting 

- — ^ ^ ^ # + ^ + ---V (35) 



dkii, \p ■ p dx dki 

The ellipses denote a term proportional to p^ , which does not contribute to the result. 
We also omit terms of order k\, which cannot contribute in the limit /c_l — >■ 0. As 
written above, the derivative acts on the hard-scattering amplitude in the momentum- 
space representation. Using an integration by parts, the derivative with respect to x can 
be made to act on the light-cone distribution amplitude. The second term, which involves 
the derivative with respect to the transverse momentum, must be evaluated before the 
collinear limit ki xp, k2 xp is taken. The light-cone projection operator of a light 
pseudoscalar meson in momentum space, including twist-3 two-particle contributions, 
then reads 



"(36) 

It is understood that, after the derivative is taken, the momenta ki and ^2 are set equal 
to xp and xp, respectively. A complete description of the pseudoscalar meson at the 
twist-3 level would also include three-particle quark-ant iquark-gluon contributions (see 
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3J] for a detailed discussion), which do not involve the large normalization factor /ip 
and thus are omitted here. (Note that the overall sign of (|36D as well as of (^if) below 
depends on the convention of ordering the quark fields in (|33|).) 

The asymptotic limit of the leading-twist distribution amplitude, valid for fi ^ oo, 
is $(a;) = 6x(l — x). For finite value of the renormalization scale, it is convenient and 
conventional to employ an expansion in Gegenbauer polynomials of the form 



^m{x,h) = 6x{l - x) 



n=l 



In numerical evaluations it will be sufficient to truncate this expansion at n 



3m and Cf''^\u) 



(37) 



2, using 



^{bu —1)- The Gegenbauer moments {n) are mul- 
tiplicatively renormalized. The scale dependence of these coefficients enters our results 
only at order a^, which is beyond the accuracy of a next-to-leading order calculation. 

The twist-3 two-particle distribution amplitudes are determined by the three-particle 
distributions via the equations of motion, except for a single term [|^. In the approxi- 
mation adopted here, where only terms proportional to hp are kept and all three-particle 
distributions are neglected, the twist-3 amplitudes must obey the equations of motion 



6 



1 — X 



%ix) 



KM 

6 



6 



(38) 



These equations enforce that we must use the asymptotic forms ^p{x) = 1 and $(t(x) = 
6x{l—x). It will be important below that ^p{x) and Ki^) do not vanish at the endpoints 
X = or 1. We finally observe that the derivative in (R3) can be substituted by 



d 



2k] 
IT 



(39) 



This is because we may expand the amplitude to first order in k± (higher powers do not 
contribute in the fc^ — limit), and use 



d 



dki 



9± ) 



(40) 



which holds after averaging k± over the azimuthal angle (denoted by {■■■))■ We use the 
definition g'^ = diag(0, —1, —1,0). The twist-3 terms in (|5B[) can now be combined into 



the projector |^| 



— — 75 



k2 ■ ki 



%{x) 



(41) 



where ki 2 are the quark and antiquark momenta defined in (0), and the factor of 
$p(x) = 1 is simply there to remind us that this is a twist-3 projection. In our anal- 
ysis below, we will quote the results of the twist-3 projections in this form, i.e., after 
eliminating (x) using the equations of motion. Expressions in terms of two functions 
$p(x) and $o-(x) are ambiguous, but reduce to the same expression upon substituting 
the asymptotic forms of the distribution amplitudes. 
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3.3 Comments on the calculation 



We now describe some technical aspects of the calculation of the various diagrams in 



more detail. The complete results for the parameters will be given in Section |3^ 



Vertex corrections 

The calculation of the four one-loop vertex diagrams (first four diagrams in Figure 
involves the trace 



tr M^^2 



(42) 



{xq + {xq + if 

where is the projector from (^), q the momentum of the meson M2, k\ (/C2) the 

momentum of the quark (antiquark) in this meson, and I the momentum of the gluon. 
We have used that M^'^^\ = ^2^*^^ = by the equations of motion, and that we can 
put ki = xq and ^2 = xq in the denominator. A Fierz transformation may be necessary 
in order to arrive at the trace (^21) . The possible F structures are therefore V — A 
(contributing to ai^... ^4^9^10), V + A (contributing to 05,7), and S + P (contributing to 

06,8)- 

li r = V ± A, only the leading-twist light-cone distribution amplitude $(x) in (RB) 



contributes under the trace. We then recover the results of [|1^, |I^ and find that there 
exist no chirally-enhanced power corrections to ai_...,5j^9_io from the vertex diagrams. If 
r = S + P, the leading-twist contribution to the trace vanishes, while the twist-3 result 
is proportional to r^^^ = 2iXM2/f^b- This gives an order correction to the coefficients 
06,8, which are multiplied by already in naive factorization. We can now exploit the 
fact that tr a^j^p = to show that the term involving the transverse- momentum derivative 
in (^) does not contribute to the trace. At this stage, the collinear limit can be taken to 
compute the kernel in the usual way. The projection on ^p{x) yields a result containing 
symmetric and antisymmetric parts under the exchange of x and (1 — x). As explained 
above, in the approximation of keeping only chirally-enhanced terms we are forced to 
assume the asymptotic form for $p(a;), so that the antisymmetric part of the kernel 
integrates to zero. The symmetric part turns out to be a scheme-dependent constant 
and is responsible for the "—6" in the expressions for ag^s in (0) below. The kernel 
resulting from the ^'„{x) projection is symmetric under x ^ (1 — x) and thus vanishes 
after integration with $^(x). 

Penguin diagrams 

We now consider the penguin contractions (fifth diagram in Figure |^), restricting our 
attention first to the twist-2 part of the projector (^). The corresponding contributions 
to the hard-scattering kernels have been given in our previous work Because there 
have been confiicting results for the penguin terms in the recent literature (see Section]^), 
we wish to clarify the origin of the discrepancies here. 

The point to note is that, depending on the structure of the {V — A) ® {V ±: A) four 
fermion operators in the effective weak Hamiltonian (without Fierz transformation!). 
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Figure 3: The two different penguin contractions. 



there exist two distinct penguin contractions with different contractions of spinor indices, 
as indicated in Figure ^. Ahhough the two diagrams are related by Fierz transformations 
in four dimensions, in dimensional regularization they give different results in the NDR 
scheme with anticommuting 75, because this scheme does not preserve the Fierz identities 
in d dimensions p^. The results for the two contractions shown in the figure involve 



left: 



2 , m; 
- In ^ 

3 /i2 



right: 



3 I 



+ 1 



(43) 



where the function G{s) is given in (^) below. The first contraction appears in matrix 
elements of the operators Q4,6, while the second one enters in matrix elements of Qi^s- 
This assumes the standard Fierz-form of the effective Hamiltonian, which is employed 
in the calculation of the Wilson coefficients in the NDR scheme. The operator 
is special, because its contribution is a pure ultraviolet effect which, by definition, is 
absorbed into the definition of the "effective" Wilson coefficient Cg^ = Csg + C5 of the 
chromomagnetic dipole operator . The discrepancies between our results and some of 
the papers discussed in Section |^ seem to arise from the fact that Q4 and Qq are treated 
like (5i,3. 

Including now the twist-3 part of the projector (|36D , we still find that the penguin 
contractions can be straightforwardly evaluated in terms of four-quark operators, before 
the actual projection is made. The on-shell conditions for the external quarks connected 
to the gluon need to be used in this step. The projection is then very simple and only 
the 75 projector multiphed by ^p{x) contributes. The resulting kernel is identical to that 
obtained with the twist-2 projection. 

The calculation of the matrix element of the chromomagnetic dipole operator (sixth 
diagram in Figure H) is more interesting in this respect. In this case it is no longer 
possible to reduce the amplitude to the usual structures involving four-quark operators. 
This is related to the fact that at the twist-3 level the fc_L momenta cannot immediately 
be put to zero. The complete twist-3 projection (|36|) has to be performed to evaluate 
the diagram, which leads to the expression 



1 



■ uMYM^'^ir^, - 7,.f )(1 + l5)uM 



(44) 



where P = Pq + k2 = Pb — ki denotes the momentum of the gluon, and the prefactor comes 
from the gluon propagator. M*^^ projects on the emission particle M2 with momentum 
q, and pb, Pq denote the momenta of the b quark and the quark in Mi, respectively. We 
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find that all four terms of the projector contribute when evaluated on this expression, 
giving the result 



ifpUqiPq) 



;i+75; 



+ Z^(l-75)k'^p(^) + o 



1 ^'(x) 1 ^Jx) 



X nif, ' ' \ 2 ^ ' ' 2 6 



+ - 



, ^(x) 2/iP , X / X 

;i +75)^ + ^(1- 75) 



X 



rrib 



X 6 

UbiPb) ■ 



Ub{Pb) 



(45) 



The second line is obtained after using the equations of motion 
distribution amplitudes, so that the asymptotic form 



for the twist-3 



1 is understood. We 



observe that the factor 1/x from the gluon propagator is cancelled in the twist-3 term, 
and so the convolution integral has no endpoint divergence. 

One can see from (^5]) that the matrix element of the chromomagnetic operator is 



obtained incorrectly at the twist-3 level if the incomplete projector containing only $p(x) 
is used. Note that, in general, this leads to gauge- dependent results, since the equations 
of motion are not respected. 



Hard spectator interaction 

The calculation of the two diagrams in the third row of Figure ^ leads to the same trace 
as in (^2|) . It therefore follows that, for T = S + P ^ the projection on M2 can result at 



most in a constant multiplying the distribution amplitude However, the constant 

obtained for the vertex diagrams resulted entirely from a term of order e = {2 — d/2) in 
the trace multiplying the ultraviolet-divergent loop diagram. Since the hard spectator 
contributions result from tree diagrams, the trace can be evaluated in four dimensions, 
and this constant is absent. We thus conclude that there is no hard spectator correction 
to the parameters 05,8 at order a^. 

If r = y ± A, only the twist-2 distribution amplitude contributes for the emission 
particle M2, but all four terms in the projector for Mi contribute to the result. We also 



find that both terms in the 5-meson projection (as given in |jT5[) contribute, but one 
of the two 5-meson light-cone distribution amplitudes drops out after implementing the 
equations of motion (0). Contrary to the other corrections, the kernels T}^ resulting 
from hard spectator interactions have logarithmic endpoint singularities at twist-3 level. 
They arise from integrals of the form J^dy/y, where y is the momentum of the antiquark 
in the meson that picks up the spectator antiquark from the B meson. These endpoint 
singularities prevent a reliable perturbative calculation of the chirally-enhanced power 
corrections to the hard spectator interactions. (The endpoint singularities are missed if 
the incomplete light-meson projector with only the 75 projection at twist-3 is employed.) 



3.4 Results for the parameters 

After these preliminaries, we now present the results for the coefficients Oj obtained at 
next-to- leading order in ag, and including the complete set of chirally-enhanced power 
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corrections to the heavy-quark hmit. We focus on the case of i? — ^ ttK decays, but similar 
results (with obvious substitutions) hold for all other B decays into two flavour-nonsinglet 
pseudoscalar mesons. For later convenience, every coefficient ailnK) is split into two 
terms: ai{7iK) = aii{7[K) + aiji{7TK). The first term contains the naive factorization 
contribution and the sum of vertex and penguin corrections (the form-factor terms in 
the factorization formula (^)), while the second one arises from the hard spectator 
interactions (the hard-scattering term in the factorization formula). Weak annihilation 
effects are not included here; they will be discussed separately in Section The 
calculation of the kernels described above results in 
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where Ci = Ci{fi), as = as{fi), Cp = {N^ — l)/{2Nc), and Nc = 3. The quantities 
V^^\ -f^MaMi) Pk2i Pk.3} Pk^i ^^^1 -P^'s^ hadrouic parameters that contain all 
nonperturbative dynamics. These quantities consist of convolutions of hard-scattering 
kernels with meson distribution amplitudes. Specifically, the terms V^^ result from the 
vertex corrections (first four diagrams in Figure H), ^'^d P^-g {P^^ and P^^) 
arise from QCD (electroweak) penguin contractions and the contributions from the dipole 
operators (fifth and sixth diagrams in Figure 0), and H^M2Mi due hard gluon exchange 
involving the spectator quark in the B meson (last two diagrams in Figure For the 
penguin terms, the subscript 2 or 3 indicates the twist of the corresponding projection. 
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In the numerical evaluation of these expressions we consistently drop higher-order 
terms in the products of the Wilson coefficients with the next-to-leading order correc- 
tions. Also, in the computation of the 0{a) corrections we confine ourselves to the 
penguin contractions of the current-current operators Qi and Q2 and to the contribu- 
tion of the electromagnetic dipole operator, which have the largest Wilson coefficients. 
These contributions, which are contained in the quantities and suffice to 

cancel the renormalization-scheme dependence of the electroweak penguin coefficients 
C7,...,io at next-to-leading order. Additional 0{a) corrections proportional to the Wilson 
coefficients C^^,,,^ of the QCD penguin operators exist, but because these coefficients are 
very small their effects can be safely neglected. 



Vertex and penguin contributions 

We now collect the relevant formulae needed for the calculation of the coefficients a^j. 
The vertex corrections are given by (M = vr, K) 



12 In 


nib 
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— X 



\nx — in 
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T / \ 1 2 ^ In X , \ 1 / \ 

2 Li2(x) — In a; H (3 + 2i7r) mx — [x ^ 1 — x) 

1 — X 



(47) 



where Li2(x) is the dilogarithm. The constants 18 and 6 are specific to the NDR scheme. 
The function g{x) can be obtained from the corresponding function relevant to, e.g., 

—>■ D^K^ decays |jl5| by taking the limit irtc —>■ 0. (Recently, a partial two-loop 
result for the vertex corrections has been obtained in , where the terms of order Poal 
were calculated analytically. To be consistent with the next-to-leading order analysis of 
this paper, we will not make use of this result.) 

If the leading-twist light-cone distribution amplitudes are expanded in Gegenbauer 
polynomials as shown in (^7^, the relevant convolution integral can be evaluated analyt- 
ically, giving 



dxg{x) ^m{x) 



?>m 

2 



Y-3.vr)«f-|af + .... (48) 



The integral with g{x) ^ g{l — x) is obtained by changing the sign of the odd Gegen- 
bauer coefficients. Since the pion distribution amplitude ^t^{x) is symmetric under the 
exchange x (1 — x), it follows that = + 12. 
Next, the penguin contributions are 
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+ (C4 + a 
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-2C: 



eff 



1 c/x 



1 — X 
i,2 -{Ci+N,C2 



pP.EW 



3ln- + --G^(s,) 



1 dx 
1 — X 



$i^(x) 



(49) 



where nj = 5 is the number of hght quark flavours, and Su = 0, Sc = {rric/mbY are mass 
ratios involved in the evaluation of the penguin diagrams. Small electroweak corrections 
from Cy, . . . , Cio are consistently neglected in 2 within our approximations. In princi- 
ple, in also contributions from C3, . . . , Cg appear. Their impact is extremely small 
numerically and we drop them for simplicity. Similar comments apply to (54) below. 
The function Gk{,s) is given by 
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Its expansion in terms of Gegenbauer moments reads 
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The contribution of the dipole operators in (^9]) involve the integral 



dx 
1 — X 



3(1 + < + «f 



+ ... 



(53) 



The twist-3 terms from the penguin diagrams are related to the twist-2 terms by the 



simple replacement ^k{x) $ 



K, 



X] 



1. For the terms proportional to C^^ and 



however, the twist-3 projection yields an additional factor of (1 
denominator in (0). We therefore find 



X] 
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which cancels the 
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with 



Gk{s) 



dxG{s - ie, 1 - x) $f (x) . 



(55) 



Inserting $^(a;) = 1, this integral is evaluated to 
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(56) 



The typical parton off-shellness in the loop diagrams contributing to the vertex and 
penguin contributions to the hard-scattering kernels is of order m;,, and hence it is appro- 
priate to choose a value /i ~ for the renormalization scale in the Wilson coefficients 
Ci and in the kernels Tj when evaluating the quantities Oj^. The //-dependent terms in 
expressions (0), (^) and ( |5^ ) cancel the renormalization-scale dependence of the Wil- 
son coefficients Cj(/i) at next-to- leading order. Similarly, the constants accompanying 
the various logarithms in these expressions are renormalization-scheme dependent and 
combine with scheme- dependent constants in the expressions for the Wilson coefficients 
to give renormalizat ion-group invariant results. 

The coefficients Ojj contain strong-interaction phases via the imaginary parts of the 
functions g{x) and G{s,x). At next-to-leading order and to leading power in KQCj^/mi), 
these phases yield the asymptotic contributions to the final-state rescattering phases of 
the B — > TiK decay amplitudes. The presence of a strong-interaction phase in the pen- 
guin function G{s^ x) is well known and commonly referred to as the Bander-Silverman- 
Soni mechanism It was included in many phenomenological investigations of non- 
leptonic B decays; however, there was always an argument as to how to choose the gluon 
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momentum, k"^, in the fourth diagram in Figure |^. In our approach there is no ambigu- 
ity to this choice, because the distribution of /c| = (1 — x)ml is determined by the kaon 
distribution amphtudes. The imaginary part of the function g{x) is another source of 
rescattering phases, which arises from hard gluon exchanges between the two outgoing 
mesons [11, |T^. The appearance of this phase is a new element of the QCD factorization 
approach. 



Hard-scattering contributions 



The hard spectator interactions shown in the last two diagrams in Figure ^ give leading- 
twist and chirally-enhanced twist-3 contributions to the kernels T/^. We include these 
hard-scattering contributions as parts of the coefficients Oj, although they are not related 
to factorized matrix elements in the usual sense. Only the twist-2 terms are dominated 
by hard gluon exchange and thus calculable. Nevertheless, for consistency with our 
treatment of the penguin coefficients we should also include the chirally-enhanced terms 
of sub leading twist, which however have logarithmic endpoint singularities. 

At tree level, the integrals over meson distribution amplitudes for the hard specta- 
tor contributions factorize. We find (the momentum fractions x and y are defined in 
Section O) 
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where we have defined the moments 



Jo 4 As 



dxx^^M{x) = (x")m • 



(57) 



(58) 



The quantity parameterizes our ignorance about the 5-meson distribution amplitude 
|T^. Not much is known about this parameter except for the upper bound SA^ < 4A 
p8| , where A = — mf, is a scheme- dependent parameter. In practice, this just means 
that \b is expected to be less than 600 MeV or so. The ratios 2fiK/mb and 



25 



multiplying the twist-3 terms coincide with the parameters and introduced in (^) 
and (§), respectively. The twist-3 contribution involves the logarithmically divergent 
integral (M = tt or K) 



■ff- 7^*;%)= T^- (59) 



lo I — y Jo 1 — y 

As previously we have to use the asymptotic form for ^^^{y), so no distinction between 
pion and kaon is necessary. (Consequently, the superscript "M" will often be dropped 
from now on.) The divergence results from the region where the spectator quark in 
the B meson enters the light final-state meson at the "lower vertex" in Figure |l] as 
a soft quark. The twist-3 hard-scattering kernels do not provide sufficient endpoint 
suppression to render this contribution subleading. In practice, the singularity will be 
smoothed out by soft physics related to the intrinsic transverse momentum and off- 
shellness of the partons, which unfortunately does not admit a perturbative treatment 
(see also the discussion in Section H). In particular, the resulting contribution may be 
complex due to soft rescattering in higher orders. For the purpose of power counting, 
we note that the effect of transverse momentum and off-shellness would be to modify 
{I — y) — > {I — y) + e with e = 0(Aqcd/''^6) in the denominator in (pU]). We thus expect 
that Xff ~ ln(mb/AQCD), however, with a potentially complex coefficient. 

Considering the off-shellness of the gluon in the last two diagrams in Figure ^ it 
is natural to associate a scale fXh ~ (Aqcd '^fc)^^^, rather than /i ~ nib, with the hard- 
scattering contributions of leading power. Hence, we set as = Cisif^h) (and also evaluate 
the Wilson coefficients at the scale fih) for the twist-2 contributions to the quantities 
Oiji in (|46|). Specifically, we use fih = V^h 1^ with A/^ = 0.5 GeV for the scale in the 
hard-scattering diagrams. However, because of the endpoint divergence the contributions 
proportional to Xff must be considered as nonperturbative effects dominated by small- 
momentum interactions. For this reason, we should more properly write 



J^Kn 



I 1\ / 1\ I ^sifJ's) ^yifJ's) I 



(60) 



and similarly for the other two quantities in (|^. Here /ig may be a soft scale. In other 
words, in principle there is no reason to expect that the twist-3 contributions should be 
governed by a perturbative coupling constant. However, it turns out that the product 
Cis{l^s) ^'xif^s) is almost renormalization-group invariant (it scales only as [as{(J?jY^'^^ with 
Nc = 3 and four flavours), and therefore we may evaluate it at the scale /ih, so that the 
ratio of running couplings in (^) equals 1. 

The expressions for the hard-scattering contributions in (^) contain many poorly 
known parameters. Besides the divergent quantity Xh and the wave-function parameter 
Xb, they depend on the 5- meson decay constant and heavy-to-light form factors. For- 
tunately, it turns out that ratios of the different hard-scattering contribution have very 
small uncertainties. Using the symmetry of the pion distribution amplitude, we flnd that 

^'ttK = ^ttK , Hktt — RttK Ht,k , (61) 
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(a) (b) (c) (d) 

Figure 4: Annihilation diagrams. 

where Rt^k = ^k-k/^t^k is the ratio of factorized matrix elements defined in ([T5|) . As ar- 
gued in Section |3.2| , in the approximation where only chirally-enhanced power corrections 
are included one is forced to employ the asymptotic forms of the twist-3 distribution am- 
plitudes ^p{x) and ^a{,x). It is then not necessary to keep other nonasymptotic effects 
in the twist-3 terms. That is, we are free to replace the moments of twist-2 amplitudes 
multiplying Xh by their asymptotic values. In this approximation, the twist-3 contribu- 
tions in (^rp reduce to a universal, multiplicative correction of the twist-2 term, and we 
obtain the approximate form of the second relation in (0). In other words, up to small 
SU(3) violations the main uncertainties in the description of the hard-scattering terms 
combine into a single poorly-known quantity H^^k- 



Finally, note that in (|57D we have not assumed any symmetry properties of the twist- 
2 light-meson distribution amplitudes. Therefore, with obvious substitutions our results 
can be applied directly to other decays, such as 5 — > tttt and Bs —>■ K^K~ . 

3.5 Weak annihilation contributions hi 

Weak annihilation contributions to charmless hadronic B decays are power suppressed 
in the heavy-quark limit and hence do not appear in the factorization formula p5|). 



Nevertheless, as emphasized in |]T^, these contributions may be numerically important 
for realistic 5-meson decays. Besides their power suppression, weak annihilation ef- 
fects differ from the hard spectator interactions discussed earlier in that they exhibit 
endpoint singularities even at twist-2 order in the light-cone expansion for the final- 
state mesons, and therefore cannot be computed self-consistently in the context of a 
hard-scattering approach. In the following discussion, we will ignore the soft endpoint 
divergences and derive results for the annihilation contributions in terms of convolutions 
of "hard-scattering" kernels with light-cone distribution amplitudes, including again the 
chirally-enhanced twist-3 projections. Despite the fact that such a treatment is not en- 
tirely self-consistent, it is nevertheless useful to estimate the importance of annihilation 
for particular final states. 

At leading order in a^, the annihilation kernels follow from the diagrams shown in 
Figure |^. They result in a further contribution to the hard-scattering term in the fac- 



torization formula, in addition to the hard spectator scattering discussed in Section ^ 
It will be convenient for future applications to keep the discussion of annihilation contri- 
butions general. We thus consider a generic 6-quark decay and use the convention that 
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M2 contains a quark from the weak decay vertex with longitudinal momentum fraction 
X, and Ml contains an antiquark from the weak vertex with momentum fraction y. The 
four-quark operators in the effective weak Hamiltonian are Fierz-transformed into the 
form (^ife)ri(^293)r2! such that the quarks in the first bracket refer to the constituents 
of the B meson. If the colour indices of this bracket are of the form qubi, diagrams (c) 
and (d) do not contribute, while diagrams (a) and (b) give rise to a colour factor Cp/Nc- 
If the colour indices are of the form qubj, then the colour factor of all four diagrams is 
Cp/N^- The projections onto the light-cone distribution amplitudes are done in the same 
manner as for the hard spectator scattering described in Section p.2| . At leading power 
(and assuming that x,y ^ ^) the integration over the 5-meson distribution amplitude is 
trivial and yields the B-meson decay constant, since the kernels are ^ independent. The 
remainders of the diagrams can be expressed in terms of the following building blocks: 
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(62) 



Here the superscripts i and / refer to gluon emission from the initial- and final-state 
quarks, respectively. The subscript k on A]l^ refers to one of the three possible Dirac 
structures TiOTa, namely k = 1 hi {V - A) ® {V - A), k = 2 for {V - A) (^{V + A), and 
k = 3 for {—2){S — P) tSi {S + P). As always, $j\/(x) denotes the leading-twist light-cone 
distribution amplitude of a pseudoscalar meson M, and the asymptotic forms of the 
twist-3 amplitudes have been used. Note that in the limit of symmetric (under x ^ x) 
distribution amplitudes, and assuming SU(3) flavour symmetry, we have A\ = A^ and 
A3 = 0. In this approximation the annihilation contributions can be parameterized by 
only two quantities {A\ and ^3). For an estimate of these annihilation contributions we 
use the asymptotic form of the leading-twist distribution amplitudes to obtain 



A\ 



Ai 



18 Xa - 4 + 



TT 



+ 2^X1 



127ra, {2X\ 



(63) 



where Xa = Jody/y parameterizes the divergent endpoint integrals. Similar to the 
case of the twist-3 hard-scattering contributions parameterized by Xh, we will treat the 
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quantity Xa as a phenomenological parameter. Clearly, taking the same value of Xa for 
all annihilation terms is a crude model. We shall see below, however, that the parameter 
aI (contributing to penguin annihilation topologies) gives the dominant contribution. 
Therefore, our treatment effectively amounts to defining a model for this particular 
parameter. 

To complete the calculation we need to account for the flavour structure of the various 
operators. It is convenient to introduce the compact notation 

(M1M2 I Ji XJ2\B,)= icfB, fM, fu, , (64) 

where the constant c takes into account factors of (—1) or appearing in the quark 

wave-functions of some of the mesons, and c 7^ only if the flavours of the "currents" ji 
and j2 match those of the mesons Mi and M2, respectively. It is apparent from Figure ^ 
that the products ji x j2 have the flavour structure 

< = E X {qiq') , (65) 

q' 

where the sum over q' = u,d,s arises from the g q'q' vertex in the annihilation dia- 
grams. Effectively, the flavour structure cr^^ "creates" a quark qi (lower index) and an 
antiquark q2 (upper index), together with a flavour- singlet q'q' pair. It is then straightfor- 
ward to find the set of meson final states to which a given flavour structure contributes. 
All a operators contributing to charged B~ decays have the structure or a". The cor- 
responding two-particle final states with light (flavour- nonsinglet) pseudoscalar mesons 
are 

a^: 71-71^, K-K^; : 7r^K~,n~K". (66) 

Operators contributing to neutral Br decays (with r = d,s) have the structure a^, a^, 
0"^, or one of the "penguin structures" tr(cr) = J2q o'q and ti^Q a) = J2q ^q<^q- Here Q 
is the charge operator for the quarks, and the sum (trace) over q is inherited from the 
QCD and electroweak penguin operators in the effective weak Hamiltonian. The final 
states to which these operators contribute are 

a^ TT+TT-, ttV, a'^: vr+i^-, 7r°ir° ; 

< : 7r°7r°, 7r-7r+, R-R^ ; 

tr((T), tr(Q a) : 7r°7r°, 7r"7r+, tt+tt", R-R+, R+R-, R^R^, R^R^ . (67) 

It follows from the above discussion that the weak annihilation contributions to the 
decay amplitudes can be parameterized in terms of the coefficients 

bi = ^ C,A\ , 63 = ^ [CsAl + C,{Al + 4) + N,C,4 
&2 = ^ C2A\ , 64 = ^ [C,A\ + CeAi 
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T EW _ 

c 

LEW _ 



C,A\ + Cr{Al + 4) + N,Cs4 



C,oA\ + CsAl . (68) 



They correspond to current-current annihilation (61,62), penguin annihilation (63,64), 
and electroweak penguin annihilation (6f^,6f^), where within each pair the two co- 
efficients correspond to different ffavour structures. The quantities 6j depend on the 
final-state mesons through the light-cone distribution amplitudes entering the expres- 
sions for A^''^ and thus should be written as 6j(MiM2). We suppress this notation when 
confusion cannot arise. As for the hard spectator terms, we will evaluate the various 
quantities in (|68|) at the scale fih = V^h yU- 

The effective weak Hamiltonian for i?-meson decays contains a strangeness-conserving 
part (TCas=o) and a strangeness-changing part {'Has=i)- Using the above definitions, the 
annihilation contributions to the matrix elements of 'Has=i can be written as 

{M,M2\nAs=i\B) = ^J2 K {M,M2\T;'^-\B) , (69) 



p= 



where Ap = VpbV*^, and 

^p"" = S^p {5rn 61 < + 5ru 63 <) + 63 < + Sm 64 tr((T) 

+ I &r er al + ^ 5rn tr(g o) . (70) 

The index r refers to the flavour of the spectator quark inside the B meson in (^), and 
6rn = Srd + ^rs squals 1 for neutral B mesons and for B^ . The matrix elements of 
TCa.s=o take an identical form, except that Ap is replaced with Ap = VpbV*^ in this case, 
and cr^ must be replaced with cx^. 

It is now straightforward to derive the annihilation contribution to a particular fi- 
nal state in terms of the coefficients hi{MiM2) . The expressions for the decay modes 
discussed in this paper have been given earlier in (|^) and (13). We will later use the 



approximation (|63|) for the quantities A^^^ to estimate the annihilation coefficients 6j 
numerically. We should recall, however, that the annihilation kernels have been derived 
under the assumption of hard scattering. Specifically, we have neglected the momentum 
fraction ^ of the spectator quark in the B meson compared to x, x, y in deriving the 
kernels. This is the reason why the results for A^^^ turned out to be independent of the 
form of the 5-meson distribution amplitude. In the endpoint regions, one or two of the 
variables x, x, y, y can be of order ^, invalidating this approximation. Therefore, our 
numerical results for the weak annihilation contributions presented in the next section 
must be considered as model-dependent estimates. 
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Table 2: Summary of theoretical input parameters. 



QCD Scale and Running Quark Masses 








m,(2GeV) 


(m„±md)(2GeV) 


225 MeV 


4.2 GeV 


(1.3 ±0.2) GeV 


(110 ±25) MeV 


(9.1 ±2.1) MeV 


Parameters Related to Hadronic Matrix Elements 


u 


h 


fB 


Fo^-(O) 


R-kK 


131 MeV 


160 MeV 


(180 ± 40) MeV 


0.28 ±0.05 


0.9 ±0.1 


Parameters of Distribution Amplitudes 












0.3 ±0.3 


0.1 ±0.3 





0.1 ±0.3 


(350 ± 150) MeV 



4 Numerical analysis of amplitude parameters 

In this section we summarize the numerical values of the parameters a, and hi entering 
the B ttK decay amplitudes and perform detailed estimates of various sources of 
theoretical uncertainties. Other decays such as i? — >• tttt will be discussed later. 

The theoretical input parameters used in our analysis, together with their respective 
ranges of uncertainty, are summarized in Table The quark masses are running masses 
in the MS scheme. Note that the value of the charm-quark mass is given at /x = m^. 
The ratio Sc = {rric/mhY needed for the calculation of the penguin contributions is scale 
independent. The values of the light quark masses are such that = r^. We hold 
{tUu + rndl/ms fixed and use as an input parameter. (This implies that in our error 
estimation procedure iPTrn/T-^nl indirectly depends on m^.) The value of the QCD scale 
parameter corresponds to as{Mz) = 0.118 in the MS scheme. The values for the 5-meson 
decay constant Jb, the semileptonic form factor F^^'^{0), and the hadronic parameter 



RjrK in (111) are consistent with recent determinations of these quantities using light-cone 



QCD sum rules form-factor models (see, e.g., [HI]), and lattice gauge theory (see, 

e.g., [0). The last row in the table contains our values for the Gegenbauer moments 
of the pion and kaon light-cone distribution amplitudes, and for the B-meson wave- 
function parameter As- The Gegenbauer moments for the light mesons are adopted with 
a conservative error estimate that encompasses most of the parameter ranges obtained 
from phenomenological or QCD sum rule determinations of these quantities. The value 
of As is an educated guess guided by the model determinations A^ ~ ~ 300 MeV 
1 3^ and Xb = (380 ± 120) MeV |3^. For comparison, using the models of Jl^ we obtain 



Xb = (410 ± 170) MeV. 
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Table 3: Next-to- leading order results for the coefficients ai^i{7rK) for three 
different choices of the renormalization scale. Numbers in parentheses show the 
maximal change in the last digit (s) under variation of the Gegenbauer moments 
of the light-cone distribution amplitudes; if present, numbers in square brackets 
show the change under variation of the charm-quark mass. 





Real Part 


Imaginary Part 




77?,;, / 2 




2?77,fc 


777,;, / 2 


' ' ^0 


2771^ 


0-1,1 


1 073f8) 


1.054(4) 


1 037f2) 


048fll) 


026f6) 


015f3) 


0-2,1 


-0 039f4) 


005f3) 


045f2) 

KJ • KJ J- '-.J \ / 


-0.113 


-0.084 


-0.066 


03,1 


0.008 


0.006 


0.004 


0.004 


0.002 


0.001 


""4,1 


0.031(4) 


0.029(3) 


0.027(2) 


0.023(0) 


0.017 


0.014 


—all 


0.036(9) [2] 


0.033(6)[1] 


0.030(4) [1] 


0.005(3) [4] 


0.004(3) [3] 


0.004(2) [2] 


-^5,1 


0.011 


0.007 


0.004 


0.005 


0.003 


0.001 




0.052 


0.052 


0.052 


0.017 


0.018 


0.019 


K c 


0.056 


0.056 


0.056 


0.005[3] 


0.007[3] 


0.008[3] 




0.007 


0.011 


0.025 


0.004 


0.002 


0.001 


r^alJa 


0.090 


0.077 


0.059 


-0.001 


-0.009 


-0.020 


r^alJa 


0.090 


0.075 


0.055 


-0.000 


-0.005[1] 


-0.010[3] 


-00,^1 a 


1.258(1) 


1.222(1) 


1.181 


0.040 


0.022 


0.012 


<o,i/a 


0.062(27) 


0.020(20) 


-0.025(16) 


0.168(39) 


0.116(29) 


0.084(22) 


o\ya 


0.062(27) 


0.018(21)[1] 


-0.028(17)[1] 


0.168(39) 


0.121(30)[1] 0.093(25)[2] 



4.1 Vertex and penguin contributions 

We begin the discussion with the vertex and penguin contributions, i.e., the terms 
ai^i{TTK) in (^). In this case all convolution integrals are finite, even for the power- 
suppressed coefficients r^aQ^iijiK) and r^a^^iij^K). Table |^ contains the values of the 
various coefficients for three different values of the renormalization scale, and with the- 
oretical uncertainties due to the variation of the pion and kaon light-cone distribution 
amplitudes and the charm-quark mass, as specified in Table |^. We vary the Gegen- 
bauer moments independently and quote the maximal variation in the table. The main 
characteristics of theoretical uncertainties are as follows: 

Renormalization- scale dependence: The residual scale dependence at next-to-leading or- 
der depends on the size of the leading-order Wilson coefficients and the magnitude of 
the Wilson coefficient that multiplies the next-to-leading order correction. In general, 
we find a significant reduction of scale dependence compared to the Oj parameters ob- 
tained at leading order (corresponding to naive factorization). Only for the parameters 
02 and aio a sizeable scale dependence remains at next-to-leading order, even though it 
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is reduced by about a factor of 2 relative to the leading order. In general, the imaginary 
parts of the coefficients a^j, which occur first at order ct^, have a larger scale dependence 
than the real parts. 

Light-cone distribution amplitudes: The explicit expressions in Section |3^ show that the 
second Gegenbauer moments af^ and enter the results for the vertex and penguin 
contributions typically with small coefficients. Therefore, the main uncertainty comes 
from the first moment of the kaon distribution amplitude, which affects only ai, 
and aiQ. The real part of oiq is particularly uncertain. There is also some uncertainty 
in 02, because the dependence on the second moment of the pion light-cone distribution 
amplitude is amplified by the large Wilson coefficient Ci. From Table |], we conclude 
that the dependence on the distribution amplitudes is almost always smaller than the 
scale dependence. Since, for consistency, we must use the asymptotic twist-3 distribution 
amplitudes, the coefficients and Og show no dependence on the Gegenbauer moments. 
This is clearly an approximation, which is valid only if the chirally-enhanced power 
corrections dominate over the remaining power corrections. (This is a questionable 
approximation, because the quark-antiquark-gluon distribution amplitude impacts on 
with a large numerical coefficient @]-) ^ consequence, we cannot control 



SU(3)-breaking effects at twist-3 order. In practice, we expect such SU(3) violations to 
have a similar (hence, small) effect as those at leading twist. 

Charm and strange-quark masses: The value of the charm-quark mass affects the penguin 
contributions with a charm-quark loop. This leads to a significant uncertainty in the 
imaginary parts of a1 and Og. The real parts of the coefficients are much less affected. 
Note that the chiral enhancement factor r^ multiplying the coefficients og and as in 
Table |^ is inversely proportional to the strange-quark mass. The ±25 MeV uncertainty in 
the value of leads to a sizeable uncertainty in the values of the products r^a^^s which, 
except for the imaginary parts of r^a^ §, is a much larger effect than the dependence on 
the charm-quark mass. 

Unknown power corrections: A source of theoretical uncertainty that is difficult to es- 
timate arises from power corrections which cannot be computed using the QCD factor- 
ization approach. Naive dimensional analysis suggests that such corrections are of order 
^QCD/nT'b ~ (10-20)%, but they could be enhanced, e.g., by large Wilson coefficients. 
This situation may potentially be realized for the parameters 02 and oiq. Recently, some 
power corrections to QCD factorization for the vrvr final state have been investigated in 
the framework of QCD sum rules and using the renormalon calculus |^6[. (Power 



corrections for final states with one heavy meson were also considered in ||15| , 0.) In 
none of these cases particularly large corrections have been identified. 



4.2 Hard spectator interactions 

We have argued in Section p.4| that the description of the hard-scattering contributions 
to the coefficients suffers from large theoretical uncertainties, which however can be 
parameterized in terms of a single (complex) quantity Ht^k defined in (|57|) . If this 
quantity is fixed, then the hard-scattering contributions a^ji can be calculated with 
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Table 4: Coefficients ai^nlnK) for three different choices of the renormahza- 
tion scale and fixed default values of the quantity H^^x- All scale-dependent 
quantities are evaluated at the scale fih = \/A^ with = 0.5 GeV. Numbers 
in parentheses show the maximal change in the last digit (s) under variation of 
the Gegenbauer moments; numbers in square brackets show the dependence 
under variation of Rttk- 





mb/2 




2mb 




1.02 GeV 


1.45 GeV 


2.05 GeV 


rrdefault 


0.92 


0.99 


1.05 




-0.087(16)[10] 


-0.061(14) [7] 


-0.045(12) [5] 


«2,II 


0.231 


0.192 


0.167 




-0.010 


-0.007 


-0.005 


'24,11 


0.004(1) 


0.003(1) 


0.002(1) 


(^5,11 


0.016 


0.010 


0.008 




-0.014 


-0.009 


-0.006 




0.112 


0.080 


0.060 


aio,ii/a 


-0.221(42) [25] 


-0.182(42)[20] 


-0.157(42)[17] 



relatively small uncertainties. In Table we show the results for these coefficients 
obtained by keeping Ht^k fixed at its central value (using central values for all input 
parameters and setting Xh = ln(mB/Aft) ^ 2.4). Because the hard-scattering terms 
arise ffist at order as, they exhibit a relatively strong scale dependence. In addition, the 
coefficients ai, and aio have some dependence on the Gegenbauer moments and on the 
value of the ratio RnK- Although the twist-3 correction is sizeable, it does not dominate 
the result for the hard spectator term. With our default value for Xh we obtain an 
enhancement of the leading twist-2 term by about 40%. 

Comparison of the results for and ji in Tables | and ^ shows that in most cases 
the hard spectator terms are of a similar magnitude as the vertex corrections. Notable 
exceptions are the coefficients 02 and aio, for which the hard-scattering contributions are 
the dominant effects. The predictions for these coefficients are correspondingly uncertain. 
On the other hand, the hard spectator contributions are very small (or absent) in the 
case of the coefficients and as. 

So far we have ignored the large overall uncertainty in the hard-scattering terms 
resulting from the uncertainty in the value of Ht^k- For an estimate of this quantity, we 
parameterize the divergent integral X^ in ( |5P| ) in the form 

Xh=(i + Qh e'""") In ^ ; Qh<1 (71) 

with an arbitrary phase ipn-, which may be caused by soft rescattering. In other words. 



34 



1 




-1 1 1 

0.5 1 1.5 2 2.5 3 

Figure 5: Ranges for the complex parameter Ht^k- The dot shows the default 
value used in obtaining the results in Table ^. 

we assign a 100% uncertainty to the "default value" Xh = ln(mA//A/j) ^ 2.4. If, in 
addition, the parameters A^, fs and Fq^^'^(O) are varied within the ranges shown in 
Table 0, the result for H^^k is confined to the interior of a region in the complex plane 
shown in Figure ^. (Variations of the Gegenbauer moments or the parameter R^^k have 
a minor effect and can be safely neglected in this plot.) The obtained values for H^^k are 
of order unity, but with an uncertainty of at least a factor 2 and a potentially significant 
strong-interaction phase (of up to about ±17° with our choice of parameters). 

4.3 Annihilation contributions 

As emphasized earlier, the results for the weak annihilation contributions derived in 
Section are based on the assumption of hard scattering, which is invalidated by the 
presence of endpoint singularities. Nevertheless, eqs. (|63|) and (^) can be employed 
as a model for the annihilation terms, which we expect to give the correct order of 
magnitude of the effects. In analogy with the previous section, we parameterize the 
divergent integral Xa in the form 

Xa= (l + f?Ae^^^)ln^; < 1 (72) 

with an arbitrary phase ipA- Table ^ shows the results for the annihilation contributions 
obtained with the default value Xa = ln(mB/A/j). They have an overall uncertainty of 
about 30% due to the error in the value of the ratio = (3.0 ± 0.9) ■ 10~^ defined in 

We observe that the default values for the annihilation contributions are rather small, 
compatible with being first-order power corrections of a canonical size. Specifically, from 
the relations for the amplitude parameters in (^) it follows that for an estimate of the 
most important annihilation effects in i? ^ ttK decays we should compare r^fea with 
(al + r^ag) (denominator of £3/2, £t, ^a), i^A{b2 + bs) with {a^ + r^a^) (numerator of Ea), 
and r^feg^ with (a^Q + r^a'^) (numerator of qc)- Similarly, from (p^ it follows that in 
B ^ nn decays we should compare r^(63 + 264) with (a^ + r'^al) (numerator of P^tt/^tttt)- 
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Table 5: Annihilation coefficients taOI for three different choices of the 
renormalization scale and fixed default values of all input parameters. 



/i 


mb/2 




2m{, 




1.02 GeV 


1.45 GeV 


2.05 GeV 


rAh 


0.025 


0.021 


0.018 


rAh 


-0.011 


-0.008 


-0.006 


rAh 


-0.008 


-0.006 


-0.005 


rAh 


-0.003 


-0.002 


-0.001 


TAhf^/a 


-0.021 


-0.018 


-0.016 


TAhf/a 


0.014 


0.010 


0.007 



In all cases, annihilation effects can be neglected in comparison with ai ~ 1. With 
the default values from Table ^ the annihilation contributions are always a moderate 
correction of less than 25% to the leading terms obtained from the QCD factorization 
formula. However, these estimates have a large uncertainty. 

Phenomenologically most relevant are the penguin annihilation effects parameterized 
by 63 and 64. They tend to increase the penguin amplitudes P in 5 — >■ i^K decays and 
PjTTr in P ^ vrvr decays, thereby reducing the values of the tree-to-penguin ratios £3/2 and 
St, and increasing the value of the penguin-to-tree ratio P,r7r/^7r7r- The fact that penguin 
annihilation graphs can significantly enhance the penguin amplitude has been noted first 
in |T^. We confirm this effect; however, with our default parameter variations we find 
a more moderate enhancement than these authors. In order to illustrate the effect and 
its dependence on the value of the quantity X^, we show in Figure ^ the combinations 
TAh and TAih + 2^4) parameterizing the penguin annihilation contributions in P — > ttK 
and B ^ nn decays. The default values for the leading penguin coefficients {a^ + r^Og) 
are shown for comparison. The regions bounded by the solid lines refer to our standard 
choice < 1 in (^) • In this case, the annihilation contribution can increase the penguin 
amplitudes by up to 30-40%. The dashed curves show the accessible parameter space in 
the more extreme case where we let Qa ^ (corresponding to a 200% uncertainty in the 
value of the divergent integral Xa)- Then the annihilation contributions can be almost 
as large as the leading penguin terms. We will see later that there are no experimental 
indications of such large annihilation effects. 

4.4 Amplitude parameters for B — ^ ttK, tttt decays 

We are now in a position to combine the results for the parameters a^j, ajji and h dis- 
cussed in the previous subsections into complete predictions for the decay amplitudes. 
We focus ffist on the amplitude parameters defined in (^Of). They are sufficient to cal- 
culate any ratio of P ^ nK decay amplitudes, such as CP asymmetries and ratios 
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Figure 6: Ranges for r^^s (left) and r^(63 + 264) (right) parameterizing, 
respectively, penguin annihilation effects in B ^ ttK and B ^ nir decays. 
Solid lines refer to qa = dashed ones to = 2. The dots show the default 
values. For comparison, the central values for the leading penguin coefficients 
(04 + r^ag) are shown by the double circles. 



of CP-averaged branching fractions. We also study the ratio PTm/T-^n in B ^ mi de- 
cays defined in (plf). A more extensive phenomenological analysis will be performed in 
Section 

It will be important to distinguish two types of theoretical uncertainties: those arising 
from the variation of input parameters to the factorization formula, and those associated 
with power corrections to factorization. Uncertainties of the first kind have a well- 
defined meaning and can, at least in principle, be reduced in a systematic way. They 
include the dependence on the renormalization scale, quark masses, moments of light- 
cone distribution amplitudes, and hadronic quantities such as and F(f~*'^(0). The 
errors in the input parameters can be reduced, e.g., by using experimental data or lattice 
calculations. The residual dependence on the renormalization scale can be reduced by 
calculating higher-order corrections to the hard-scattering kernels. Our predictions also 
depend on the value of iV^b/Kfeh however, this should not be considered a theoretical 
uncertainty. Ultimately, the goal is to use hadronic B decays to learn about such CKM 
parameters. 

Theoretical uncertainties related to power corrections to the factorization formula are 
of a different quality. Since factorization does, in general, not hold beyond leading power, 
there is no systematic formalism known that would allow us to analyze power corrections 
in a model-independent way. This is a general problem of QCD factorization theorems in 
cases where no operator product expansion can be applied. (Another familiar example 
are event-shape variables in e'^e~ annihilation to hadrons.) In the present work, we have 
identified sources of potentially large power corrections (chirally-enhanced contributions 
and weak annihilation terms) and estimated their effects. These estimates are uncertain 
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due to logarithmically divergent endpoint contributions, which indicate the dominance 
of soft gluon exchange. Without significant conceptual progress in the understanding of 
power corrections to observables that do not admit a local operator product expansion, 
it will be difficult to reduce these uncertainties in a systematic way. 

Our results for the amplitude parameters including all theoretical uncertainties of the 
first kind are shown in Table ^. They are obtained by keeping the parameters Xh and 
Xa entering the power corrections fixed at the default value Xh = Xa = \n.{mB / ^h)- All 
other input parameters are varied within the ranges shown in Table |^. Following common 
practice, we vary the renormalization scale /i between nib/ 2 and 2m},. The individual 
contributions to the error are then added in quadrature to obtain the total theoretical 
uncertainty of the first kind. The two most important contributions to the total error are 
shown in the last column of the table. The sign convention is such that the upper (lower) 
sign corresponds to increasing (decreasing) the value of an input parameter. Finally, the 
second error on the central value (if present) indicates the sensitivity to the uncertainty 
in the ratio = 0.085±0.017, which we assume to be 20%. This is not a hadronic 

uncertainty and therefore should not be combined with the first error. The amplitude 
parameters are either proportional (or inversely proportional) to |14b/K6| or independent 
of this parameter, so that this error can easily be readjusted if needed. 

In almost all cases the theoretical uncertainty is dominated by a single source. With 
the exception of the strong-interaction phase ujc, our next-to- leading order results for 
the amplitude parameters are very stable under variation of the renormalization scale. 
The uncertainty in the values of the tree-to-penguin ratios £3/2, ^r, and |P^^/T^^| is 
dominated by the error on the strange-quark mass, whereas the corresponding strong- 
interaction phases 0, 0^, and arg(P^7r/^7r7r)7 as well as the phase 0^, are most sensitive 
to the error on the charm-quark mass. The largest uncertainty in the electroweak- 
penguin parameter q comes from the SU(3) violations parameterized by the amplitude 
ratio Rt^k in ([T3|). The uncertainty in the value of the wave-function parameter \b is 
the dominant source of uncertainty for the electroweak-penguin parameter qc- Note 
that the Gegenbauer moments of the pion and kaon light-cone distribution amplitudes 
are never the dominant contribution to the error. This shows that the precise shapes 
of these amplitudes are of minor importance for phenomenological applications of QCD 
factorization. 

In order to study the theoretical uncertainties of the second kind, related to our es- 
timate of potentially large power corrections to the factorization formula, we show in 
Figure |^ results for amplitude parameters obtained using central values for all input 
parameters, but varying the complex quantities Xh and Xa according to the param- 
eterizations ([7l|) and (|72D . The dots with error bars show the central results for the 
amplitude parameters as given in Table ^ (for fixed |\4fe/K:fe| = 0.085). The dashed 
curves (visible only for £3/2 and qc) bound the parameter space obtained by variation 
of Xh-i whereas the solid curves (present for all parameters except q) define the region 
obtained by variation of Xa- It is evident that the chirally-enhanced twist-3 corrections 
to the hard spectator interactions do not lead to a dominant uncertainty. Their effect is 
always smaller than the uncertainty due to parameter variations of the first kind. The 
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Table 6: Predictions for the amplitude parameters including all theoretical 
uncertainties of the first kind (sec text), but using default values for the power 
corrections to QCD factorization. The first error on the central value is the 
sum of all theoretical uncertainties added in quadrature. The second error (if 
present) shows the dependence on \ Vub/Vch\- The last column indicates the two 
most important contributions to the theoretical uncertainty. For each quan- 
tity, the second line shows the result without weak annihilation contributions 
(except for q and u, which do not receive annihilation terms). 



Parameter 


Central Value 


Dominant Errors 


£3/2 (%) 
(j) (deg) 


23.9 ±4.5 ±4.8 
25.7 ±4.8 ±5.1 
-9.6 ±3.8 
-10.2 ±4.1 


±3.5 (m«) ±1.4 (/i) 
±3.6 (m,) ±1.6 («f) 
±3.5 (mj ±1.4 (af) 
±3.7 (mj ±1.5 (af) 


St (%) 
0T (deg) 


20.6 ±3.5 ±4.1 
22.0 ±3.6 ±4.4 
-5.7 ±4.4 
-6.2 ±4.6 


±3.2 (m,) ±0.9 (/i) 
±3.3 (m,) ±0.8 (af ) 
±3.5 (mj ^2.3 (/x) 
±3.7 (mj T2.2(/i) 


Sa (%) 

0a (deg) 


2.0 ±0.1 ±0.4 
1.9 ±0.1 ±0.4 
13.6 ±4.4 
16.6 ±5.2 


±0.1 (mc) t0.1(/x) 
±0.1 (mj — 
±3.7 (me) ±1.7 (af) 
±3.9 (me) t2.8 (/i) 


q{%) 
uj (deg) 


58.8 ± 6.7 T 11.8 
-2.5 ±2.8 


±6.4 {R^k) ±1-3 (/i) 
±1.9 (/x) Tl.8(af) 


qc (%) 
uc (deg) 


8.3 ± 4.5 T 1-7 
8.9 ± 4.9 T 1.8 
-60.2 ±49.5 
-54.2 ±44.2 


T2.7(Ab) ±2.3 («f) 
T3.1(Ab) ±2.3 (af) 
±31.7 (/x) ^27.9 (Ab) 
±29.5 ill) ^24.1 (Ab) 


l-fTTTr/^TTTrl (9o) 

arg(P^^/T^^) 

(cleg) 


28.5 ±5.1^5.7 
25.9 ±4.3^5.2 
8.2 ±3.8 

9.0 ± 1.1 


:f4.6 (m,) ^1.8 
:F4.1(m,) =F0.8(At) 
T3.3 (me) ±2.0 (//) 

T3.6(///e) ±1.8(//) 
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Figure 7: Results for the amplitude parameters. Dots with error bars show 
the default values and uncertainties of the first kind, while the regions bounded 
by the curves show the variation of the central values under variation of Xh 
with 0H = ^ (dashed curves), and variation of Xa with Qa = ^ (solid curves). 



uncertainty in the description of weak annihilation is significant for the tree-to-penguin 
ratios £3/2, St, Pmr/Tn-K and their phases. The resulting variation is typically as large as 
the total theoretical uncertainty of the first kind. Although the modelling of weak anni- 
hilation thus introduces sizeable uncertainties, we stress that their impact on the values 
of the amplitude parameters is still a moderate correction (see also Table P). Moreover, 
it is important that the dominant effect is a universal contribution to the leading pen- 
guin amplitudes P and Ptttt, which leads to correlations between the various amplitude 
parameters. Ultimately, this will help to constrain the annihilation contributions using 
experimental data (see Section below). 



5 Phenomenological applications 

This section, which can be read without studying the technical details of our analysis, 
illustrates several applications of our result to the phenomenology of 5 — nK and B 
vrvr decays. The main goal will be to obtain information about the Wolfenstein parameters 
p and fj defining the unitarity triangle. After some general remarks in Section we 



40 



start by considering methods for constraining p and fj that depend on a minimal amount 
of theoretical input about QCD dynamics. The simplest such strategy is the Fleischer- 
Mannel bound P] discussed in Section f).2\ In the following Section p.3|, we review 



analysis strategies developed by Rosner and one of us fl], H @], which are based on 
CP-averaged rate measurements for the charged modes -B^ — > nK^ vrvr. These methods 
provide powerful constraints in the (p, ff) plane with minimal theoretical uncertainties. 
The third method, the determination of sin 2a from the time-dependent CP asymmetry 
in , —>■ TT+TT^ decays, is explored in Section |5]^. Towards the end of our discussion 
we will then rely more heavily on the new theoretical results obtained in the present 



work. In Section |5.5| , we give predictions for the absolute values of branching fractions 
and results for various ratios of CP-averaged branching ratios as a function of 7. We 
then perform a global fit to the experimental data on the branching ratios and extract 
the corresponding allowed region in the (p, fj) plane. In the final Section |5.6| , we show 
predictions for the direct CP asymmetries in the various B ttK and B ^ tttt decay 
modes. Because of their sensitivity to strong-interaction phases, these results have the 
largest theoretical uncertainties. 



5.1 General observations 

Several points following from the numerical analysis in Section |4.4| are worth repeating 
here, because they will have a direct impact on some of the analysis strategies mentioned 
below: 

1. The rescattering effects parameterized by Ea are very small and not much affected 
by theoretical uncertainties. Therefore, the decay amplitude for B~ — >■ n~K^ in (0) is, 
to a very good approximation, given by the pure penguin amplitude P. This has two 
important consequences: first, it is a safe approximation to neglect terms of order in 
the squared decay amplitudes; secondly, it follows that the direct CP asymmetry in the 
decays B^ 7i^K°, 

Acp(7r+ir°) = -2ea sin (j)a sin 7 + 0{el) ^ -1% x sin 7 , (73) 

is tiny and unobservable in the foreseeable future. (We define the CP asymmetries as 
the difference of the 5- meson minus S- meson decay rates divided by their sum.) An 
experimental finding of a sizeable asymmetry in this decay mode would have to be 
interpreted as a sign of physics beyond the Standard Model, or as a gross failure of 
the QCD factorization formula, indicating the presence of large, uncontrollable power 
corrections. For completeness, we note that the prediction that is small can be 
tested experimentally by measuring the CP-averaged branching ratio for the decays 
5± K^K^ [I 15]. 

2. The strong-interaction phase u = —(2.5 ± 2.8)° is accurately predicted and tiny, 
consistent with zero within errors. Also, the value of q is not affected by annihilation 
contributions. These observations confirm a theoretical argument presented in @, |^, 
which uses Fierz identities and top-quark dominance in the electroweak penguin diagrams 
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to show that ~ 0, and that q can be calculated in a model- independent way up to 
SU(3)-breaking corrections. The argument is based on the fact that in the limit of V-spin 
(s ^ u) symmetry the leading electroweak penguin contributions parameterized by q e^^ 
can be related to the current-current contributions from the operators and Q'2 using 
Fierz identities. Neglecting the small contributions from the operators Q7, Qs? and Qy^, 
one then obtains 



where the ratio of Wilson coefficients is renormalization-scheme invariant and can thus 
be evaluated at the electroweak scale. The strong-interaction phase uj vanishes in this 
approximation. When SU(3)-breaking corrections and small electromagnetic contribu- 
tions are included, the above result is rescaled by a factor Rq = (0.84 ±0.10) _ 
About half of the deviation from 1 is due to (mostly "factorizable" ) SU(3) violations. 
The important point to note is that the smallness of is a model-independent result 
that does not rely on the QCD factorization formula. It follows that terms of second 
order in u can be safely neglected. 

3. The electroweak penguin contributions to the 7i~^K~ decay amplitude, param- 
eterized by qc in (0), are small and can also be treated to first order to good approxi- 
mation. In the literature, these effects are sometimes referred to as "colour-suppressed" 
electroweak penguin contributions. 

4. The strong-interaction phases 0, 4>t, and arg(P^^/T^^) are small. We find central 
values of order 10° or less in magnitude, with an uncertainty of about a factor 2 due 
to potentially large annihilation contributions and higher-order perturbative corrections 
to the hard-scattering kernels. It follows that the cosines of these phases deviate from 
1 by only a few percent, and the direct CP asymmetries are suppressed by a factor 
|sin0j| ~ 0.1-0.3. We note, in this context, that the smallness of the phase difference 

— 4° (see Table |]) can be understood based on simple physical arguments and 
implies a strong correlation between the direct CP asymmetries in the decays 



In the following sections, we explain how these general observations can be put to work in 
different analysis strategies. We start with those strategies that avoid theoretical input 
on tree-to-penguin ratios such as £3/2 and Et- This eliminates the sensitivity to weak 
annihilation effects and hence the main uncertainty of the QCD factorization approach. 
As a result, these strategies are particularly clean from a theoretical point of view. 

Experimental data for the CP-averaged B ttK and B ^ Tin branching fractions as 
reported by several experimental groups are collected in Table |^. The last column shows 
our naive averages of these results neglecting correlations. We use the average result for 
the B^ — > 7r^7r° branching ratio in despite of the fact that the individual measurements 
of this mode have less than 3cr significance. 




(74) 



TT^K^ and B' 
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Table 7: Experimental results for the CP-averaged B — ^ ttK and B ^ tttt 
branching ratios in units of 10~^. The BaBar and Belle results are preliminary. 
Our averages ignore correlations. 



Decay Mode 


CLEO 47] 


BaBar 4| 


Belle [49| 


Average 




4.31};^ ±0.5 
5.6l|^±1.7 
< 12.7 (90% C.L.) 


4.1 ±1.0 ±0.7 
5.11?;^ ±0.8 
< 9.0 (90% C.L.) 


5.9i|t ±0.5 

7 1 +3.6+0.9 
' -J- -3.0-1.2 

< 12.6 (90% C.L.) 


4.4 ±0.9 
5.6 ± 1.5 




17.2t|^± 1.2 
11 +3.0+1.4 

ii.D_2. 7-1.3 

18.2tto± 1-6 


16.7± 1.61J:? 

1 n 0+2.1+1.0 
J-U.o_1.9_1. 2 

1 0+3.3+1.6 

-'-"•^_3.0-2.0 

9+3.1+1.1 
o-^ -2.7-1.2 


18.7l||j±1.6 

1 7 n +3.7+2.0 
J- '-u -3.3 -2.2 

13.lt|^±2.6 
14.6l^:J ±2.7 


17.2 ± 1.5 

12.1 ± 1.7 

17.2 ±2.5 

10.3 ±2.5 



5.2 Fleischer— Mannel bound 

Based on a few plausible assumptions, Fleischer and Mannel have proposed a very simple 
method for obtaining a bound on the weak phase 7 from the measurement of a single 

_ t{B+) Br(E° ^ 7r-ir+) ± Br(B° ^ Tr+jiT-) 

~ 7(5oy Br(S+ ^ TT+i^o) ± Br(5- ^ 7r-i?0) ^ ^ 

of CP-averaged branching ratios (see for more sophisticated generalizations of this 
method). The current value of this ratio obtained from Table ^ is i? = 1.06 ± 0.18. 
Neglecting the small rescattering effects parameterized by Sa and the colour-suppressed 
electroweak penguin contribution qc (observations 1 and 3), it follows from (|18|) that 

i? ~ 1 — 2eT cos cos 7 ± el > sin^ 7 . (76) 

The bound excludes a region near I7I = 90° provided that R < 1. It is valid for any 
values of Et and the strong-interaction phase (pT and is thus independent of theoretical 
assumptions about the tree-to-penguin ratio in these decays. 

Our calculations confirm the dynamical assumptions that go into the derivation of 
this bound. The neglected terms in the exact formula for R are of order EaSx or QcSt- 
Both are at the few percent level and can be safely neglected until the experimental error 
on R is much below 10%. (More accurately, we find that the bound ([76|) is violated by 
at most 1.5% and only in the region 72° < 7 < 86°.) 

5.3 Strategies based on charged modes 

The theoretical analysis of the charged decays B^ ttK profits from the fact that, with 
the exception of the strong- interaction phases and (pa, the hadronic parameters entering 
the parameterization of the corresponding decay amplitudes in (|18|) can be constrained 
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in the limit of SU(3) symmetry 0, ^, ^. Therefore, the QCD factorization formula is 
needed only to reduce the uncertainties in the estimate of SU(3)-breaking corrections. 
We have already seen above how SU(3) symmetry and Fierz identities help to calculate 
the quantity q e*"^ with small theoretical uncertainties. In addition, the decay amplitudes 
for the charged modes — > nK in ([T8|) depend on the tree-to-penguin ratio £3/2 and 
the very small rescattering parameter (and the corresponding phases (p aiid (pa)- It 
is apparent from Figure |^ that our prediction for £3/2 has a large uncertainty due to 
weak annihilation contributions as well as parameter variations. The advantage of the 
charged 5^ — * ttK modes is that the tree-to-penguin ratio can be determined, up to 
small SU(3) violations, using experimental data. Specifically, a certain combination of 
the parameters £3/2 and Ea can be measured by comparing the CP-averaged branching 



fractions for the decays B 



£3/2 



± 



7r±7r° and B^ 



'^K'^. The relation is 



^3/2 



1 + 2ea cos (pa COS 'J + 



exp ) 



where 



'Cxp 



tan 6c ^ 

Jtt 



2[Br(fi^ 



1/2 



(77) 



(78) 



is an observable, and Rth = 1 in the limit of U-spin [d ^ s) symmetry 0. In the 
theoretical analysis of irK decays it is convenient to replace the parameter £3/2 by 

£3/2. Since = O(eKM) is very small (observation 1), both quantities take very similar 
values. The experimental uncertainty in the current value e, 



exp 



0.223 ± 0.034 is still 

large. Ultimately, however, the accuracy in the determination of £3/2 is only limited 
by the theoretical uncertainty in the calculation of the SU(3)-breaking corrections to 
Rth- The QCD factorization approach helps reducing the model dependence in this 
calculation. Neglecting the tiny contributions from electroweak penguins. 



R 



th 



ai{TTK) + R^K 0.2(71"-^) 



ai(7r7r) + a2(7r7r) 



m 



K) 



(79) 



K 



The form-factor ratio can be safely set equal to 1. (Deviations are of order 
m^)/m?Q 1%.) Note that there are no annihilation contributions to -Rth- If we keep 
the parameter governing the twist-3 contributions to the hard spectator interactions 
fixed vary all other input parameters over their respective ranges of uncertainty, we find 
that Rth = 0.98 ± 0.02. Next, if we keep the input parameters fixed but vary Xh as 
shown in ([7I|), we find a variation of about ±0.01. However, since our focus here is on 
SU(3) violations we should be more conservative and let the quantities X^ and vary 
independently. This gives a larger variation of about ±0.03. Taken altogether, we obtain 

Rth = 0.98 ± 0.05 . (80) 

Combining this with the experimental value of ^exp gives £3/2 = 0.218±0.034cxp±0. 011th- 
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Figure 8: Comparison of the prediction for the tree-to-penguin ratio £3/2 in 
5^ — ^ ttK decays with the experimental value (at la and 2a) of the quantity 
£exp defined in (|7Sp. The region bounded by the solid line refers to our default 
annihilation model with qa = ^- The dashed line corresponds to qa = 2. 



Our results for the SU(3)-breaking corrections parameterized by -Rth and Rq (the 
quantity that corrects (0)) are valid up to nonfactorizable corrections that simultane- 
ously violate SU(3) symmetry and are power-suppressed in Aqcd/''^6- The potentially 
most important power corrections are included in the error estimate. The remaining 
uncertainties are of order 

^ / 1 m, — , , 

0[^-— 81 

ViV, mb y" ^ ' 

and, by naive power counting, should not amount to more than a few percent. 

Whereas our theoretical predictions for the parameter £3/2 were affected by large un- 
certainties (see Figure |^), relations ( [77| ) and ( pUD can be combined to obtain a much more 
accurate value for the related parameter £3/2 provided, of course, that the experimental 
value of ^exp has a small error. This parameter can then be used in the phenomeno- 
logical analysis of ttK decays. Moreover, the comparison of the so-determined 

value of this parameter with our prediction for the tree-to-penguin ratio £3/2 provides a 
nontrivial test of the QCD factorization approach, and ultimately could help to reduce 
the uncertainties in our modelling of weak annihilation contributions. 

To illustrate this last point, we show in Figure ^ our prediction for £3/2 e**^ (corre- 
sponding to the first plot in Figure 0) and underlay as a gray band the experimental 
value for ^exp with its la and 2a errors. (The two quantities should agree if we neglect 
the small deviation of i?th from 1, and the tiny contribution of Sa in (^).) It is pleasing 
that the central experimental value is in excellent agreement with our theoretical pre- 
diction. The "smallness" of the tree-to-penguin ratio in 5 — nK decays has often been 
interpreted as evidence for large, nonfactorizable contributions to the decay amplitudes. 
Here we find that this effect is reproduced in the QCD factorization approach without 
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any tuning of parameters. This is an important result. A deviation of our prediction 
from the experimental result could have been considered as an indication of large cor- 
rections to QCD factorization, such as enhanced weak annihilation effects not covered 
by our simple model estimates. As an example, the dashed curve in the figure shows the 
allowed region obtained by increasing the value of the parameter qa in (|72[) from 1 to 2. 
Fortunately, the data provide no evidence for the existence of such large deviations from 
our central prediction. (However, the evidence from £exp alone does also not exclude a 
potentially large annihilation contribution, if it has a large phase.) 



Model-independent bound on 7 

A key observable in the study of the weak phase 7 is the ratio of the CP-averaged 
branching ratios in the two -B^ ttK decay modes, defined as 

Br(S+ ^ TT+fsfO) + Br(fi- ^ n-R^) 
* ~ 2[Br(fi+ ^ nOK+) + Br(fi- ^ nOK+)] ' ^^^^ 

Its current value is i?* = 0.71 ± 0.14. The theoretical expression for this ratio obtained 
using the parameterization in (|T8|) is 



i?/ = l + 2£:3/2Cos0(g-cos7)+£:3/2(l-2gcos7-Fg^) 

— 2^3 /2 €a sin^7 COS (p COS 0a + ( 1 — g COS 7) sin sin 0^ 

- 2^3/2 q sin sin + 0{el, , Ea uj) . (83) 

(The parameter q is also called 5ew in the literature on this ratio.) Note that the 
rescattering effects described by as well as the terms linear in the small phase cu, are 
suppressed by a factor of £3/2 and thus reduced to the percent level. Higher-order terms 
in these small parameters can be safely neglected (observations 1 and 2). 

From a measurement of the ratio i?* alone a bound on cos 7 can be derived, which 
for i?* < 1 implies a nontrivial constraint on the Wolfenstein parameters p and f/ 0]. 
Only CP-averaged branching ratios are needed for this purpose. The idea is to allow 
the strong-interaction phases and 0^ in ( |5BD to take any value between and 27r. 
(Of course, the QCD factorization approach predicts rather restricted ranges for these 
strong-interaction phases. For the moment, however, we will not make use of these 
predictions.) This gives p, || 

< {1 + e3/2\q - cos'jiy + €3/2(63/2 + 2ea) sin^-f + 0{el,€aUJ,uj^) . (84) 

Note that there is no term linear in uj on the right-hand side. Provided R^, is significantly 
smaller than 1, the bound implies an exclusion region for cos 7, which becomes larger 
the smaller the values of i?* and £3/2 turn out to be. The effect of the rescattering 
contribution proportional to Ea on the right-hand side of the bound is numerically very 
small. 
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Figure 9: Theoretical constraints on the Wolfenstein parameters (p, ff) imphed 
by (still hypothetical) experimental upper bounds on the ratios i?* and ^exp- 
The left-hand plot refers to ^exp < 0.19, the right-hand one to £exp < 0.25. For a 
given upper bound on i?*, p and fj are restricted to lie inside the corresponding 
shaded region. The dashed circles show the allowed region implied by the 
measurement of \Vub/Vcb\ in semileptonic B decays. 



For fixed value of g, eq. (|^) excludes a region in cos 7 provided that <1. How- 
ever, we should take into account that the value of q itself depends on the Wolfenstein 
parameters p and fj. Hence, it is more appropriate to display the relation ( p4D as a 
constraint in the (p, fj) plane. We use 

p q 0.222 ±0.025 

VP + Tj VP + V 

where the numerical value q = q Rb = 0.222 ± 0.025 corresponds to iV^t/Kfel = 0.085 ± 
0.017. The experimental inputs to the bound are the measured ratios i?* and ^exp of CP- 
averaged branching ratios. The theoretical inputs are the value of q, the parameter Rth in 
the relation £3/2 = i?th^exp, and a value for the rescattering parameter Sa- The accuracy 
of our predictions for q and i?th is intrinsically limited only by effects that have a strong 
parametric suppression, as shown in (|8T|). As far clS S(i IS concerned, the bound becomes 
weaker the larger the value of Ea- In our analysis we take Ea < 0.04, corresponding to an 
upper bound that is twice as large as predicted by the QCD factorization approach. 

Figure ^ illustrates the resulting constraint in the (p, fj) plane obtained for some 
representative upper bounds on R^ and ^exp- For comparison, the dashed circles show the 
constraint arising from the measurement of the ratio \ Vub/Vcb\ in semileptonic B decays. 
(Note that the information from kaon CP violation excludes f/ < in the Standard 
Model.) It is evident that, depending on the values of i?* and ^exp, the constraint may be 
very nontrivial. If R^, < 0.7, then only values I7I > 90° are allowed, which are significantly 
larger than those favoured by the global analysis of the unitarity triangle (see ^ 



47 



52| , |53[] for some recent discussions of the standard analysis). For yet smaller values 
< 0.6, the arising constraint would become inconsistent with the global analysis. 
This would be an indication of some new flavour physics beyond the Standard Model. 
On the other hand, if i?* > 0.9 then the excluded region is too small (or absent) to be 
of phenomenological significance. The present uncertainty in the value of i?* is too large 
to tell which of these possibilities is realized. 

Determination of 7 and constraint in the (p, fj) plane 

Ultimately, the goal is of course not only to derive a bound in the (p, ff) plane, but to 
determine the Wolfenstein parameters (and thus the unitarity triangle). This requires 
obtaining information about the strong-interaction phase in (|83D , which can be achieved 
either through the measurement of a CP asymmetry or with the help of theory. A 
strategy for a model-independent determination of 7 from — * vrfC, vrvr decays has 
been suggested in [ffl]. It generalizes a method proposed by Gronau, Rosner and London 



5J] to include the effects of electroweak penguins. The approach has later been refined 
to account for rescattering contributions to the 71"^ decay amplitudes 

However, this method relies on the measurement of a direct CP asymmetry in addition 
to i?* and hence requires very high statistics. Here, we suggest an easier strategy for a 
theory-guided determination of p and f], which does not require the measurement of a CP 
asymmetry. Instead, we will exploit the prediction of the QCD factorization approach 
that the strong- interaction phase is small, i.e., 

sin = 0[a^(mfe), AqcoMb] . (86) 

This implies that the deviation of cos0 from 1 is a second-order effect in as(mb) and/or 
Aqco/'^b- More specifically, we found that ~ —(10 ± 15)°, where the uncertainty is 
dominated by the weak annihilation contribution. Taking this value literally we would 
conclude that cos0 > 0.9; however, to be conservative we allow for a larger phase such 
that COS0 > 0.8. (This prediction can be tested experimentally once the direct CP 
asymmetry in the decays B^ tt^K^ has been measured |1|, ^.) With the help of this 
result, a measurement of the ratio i?* can be used to obtain a narrow allowed region in 
the {p,fi) plane, which for fixed value of |V^b/K;fe| corresponds to a determination of 7 
that is unique up to a sign. To this end, we rewrite (^) as 

cos7 = g + — — / ^ + 0{e^,eau;,uj ) , (87) 

263/2 (cos £3/2 q) 

where 

5 = Ea sin^7cos0cos0a + (1 — gcos7) sin0sin0a + g sin sin c<j . (88) 

It is safe to take < 6 < 0.05 and treat this as an independent parameter. Then, in 
addition to q and -Rth, we must specify a range for cos0. 

In evaluating the result (|87D , we scan the theory input parameters in the ranges 
0.197 <q< 0.247, 0.93 < i?th < 1-03, <6 < 0.05, and 0.8 < cos0 < 1 (corresponding 
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Figure 10: Left: Allowed regions in the {p^ff) plane corresponding to ^exp = 
0.22 and different values of i?* as indicated. The widths of the bands reflect 
the total theoretical uncertainty, obtained by scanning the input parameters 
inside the ranges specified in the text. Right: Allowed regions at 68% and 
95% confidence level obtained from the current experimental results on the 
branching ratios. The dark band shows the theoretical uncertainty for the 
central experimental values. The light gray area is the allowed region obtained 



from the standard global fit of the unitarity triangle p3| 



to 101 < 37°). The left-hand plot in Figure TO shows the allowed regions in the (p, ff) plane 
obtained for some representative values R,, and the current central value ^exp = 0.22. We 
stress that with this method a useful constraint on the Wolfenstein parameters is obtained 
for any value of these parameters. When combined with a measurement of iVub/Kfel? this 
determines the weak phase 7 up to a sign ambiguity. Note that the theoretical accuracy 
of the method is high, especially if the ratio turns out to be close to 1, corresponding 
to a value of 7 as suggested by the global analysis of the unitarity triangle. In that 
case the theoretical uncertainty in the determination of 7 is about 10° (for fixed value of 
iKfe/KbD- Also, the resulting constraint is then very weakly dependent on the value of 
the parameter ^exp- We stress that the width of the band in the (p, ff) plane corresponding 
to i?* = 1 is narrower than the widths of the theoretical error bands corresponding to 
the "standard" constraints on the unitarity triangle derived from charmless semileptonic 
B decays, B-B mixing, and CP violation in K-K mixing. Only the measurement of 
sin 2/3 is theoretically cleaner. The uncertainty in the extraction of 7 increases as one 
considers values of i?* significantly less than 0.8 or larger than 1.2. However, such values 



would be inconsistent with the global analysis of the unitarity triangle [^3) 11) e1 
and thus provide evidence for physics beyond the Standard Model. 

With the current values of the branching ratios collected in Table |^, the method 
proposed here is at the verge of providing a useful constraint in the (p, ff) plane. This is 
shown in the right-hand plot in Figure |l^, where we indicate the resulting allowed regions 
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at 68% and 95% confidence level and compare them with the allowed region (light gray 
area) obtained from the standard global fit of the unitarity triangle. Here and below we 
use the most recent result for the standard fit obtained in [^], which includes the mea- 
surements of sin 2(3 at the i?-factories and adopts a conservative treatment of theoretical 
uncertainties that is similar in spirit to the one adopted here. In evaluating the la and 
2cr domains of the quantities R,, and ^exp? we take into account the correlation implied 
by the fact that the —>■ tt^K^ branching ratio enters both quantities. Specifically, 
we vary the branching fractions for —>■ tt^tt'^, tt^K^ and n^K^ independently such 
that ^ 1 (for Icr) o^' 4 (for 2a). In this way we find that the minimum and maximum 
values of i?* at 95% confidence level are 0.47 and 1.07, respectively. If in the future the 
upper value can be reduced, the resulting allowed region in the (p, f/) plane will no longer 
fully overlap with the standard domain. 

5.4 Determination of sin 2cx from B — > tt+tt" decays 

The methods described so far in this section provide constraints in the (p, fj) plane that, in 
essence, correspond to a determination of the magnitude of 7 = arg(\4*^). Independent 
information about the unitarity triangle can be obtained from a measurement of the 
time-dependent CP asymmetry in the decays B^, B^ tt+tt^, which is sensitive to the 
Bd-Bd mixing phase e~^*^. We define 



Br(SO(t) ^ TT+TT-) - Br(SO(t) ^ tt+tt-) 
Br(50(t) ^ TT+TT-) + Br(50(t) ^ n+n-) 

— S'tttt sin(AmB t) + Ctttt cos(Ams t) , (89) 



where 



Q 2ImA,r7r ^ 1 IAttttI , —2if5 ^ "I" ^TTTr/^TTTr (Qf\'\ 

~ 1 I I X 12 ' ^'^'^ ~ 1 I I \ 12 ' '^'^'^ ~ ^ pn I p Jt ■ ^ ' 

± ~r I ^TTTT I ~r I /Vti-tt I O "T J TTTT/ TTTT 

The coefficient Ctttt, which is a function of the weak phase 7, coincides with the direct CP 
asymmetry to be discussed later. The mixing-induced asymmetry S'tttt depends on 7 and 
f3. In fact, in the limit where P^^/T^^ is set to zero it follows that A^^ = e~^*^^+''') = e^*", 
and hence Sj^t, = sin 2a. 

When the penguin contributions to the B ^ decay amplitudes are included, the 
relation between the coefficient St,tt and sin 2a receives hadronic corrections |P, |^, plj] . 



which can be calculated using the QCD factorization approach To illustrate the 

effect, we first assume that and the weak phase /5 have been determined accu- 

rately. Then using 7 = 180° — a — /3 the expression for A^j^ in ( pOD becomes a function of 
a and our prediction for the penguin-to-tree ratio P^^/T^tt- If we further assume that the 
unitarity triangle lies in the upper half of the (p, f/) plane, then a measurement of Smr de- 
termines sin 2a with at most a two-fold discrete ambiguity. Figure |Il| shows the relation 
between the two quantities for the particular case where \Vub/Vcb\ = 0.085 and /? = 14.3°, 
corresponding to sin 2/? = 0.48 (the current world average is sin 2/3 = 0.48 ± 0.16 
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Figure 11: Relation between sin2Q; and the mixing- induced CP asymmetry 
5*^^, assuming sin 2/5 = 0.48. The dark band reflects parameter variations 
of the first kind, the light band shows the total theoretical uncertainty. The 
lower portion of the band refers to values 45° < a < 135°, the upper one to 
< a < 45° (right branch) or 135° < a < 180° - /3 (left branch). 



The dark band shows the theoretical uncertainty due to input parameter variations as 
specified in Table 0, whereas the light band indicates the total theoretical uncertainty 
including the effects of weak annihilation and twist-3 hard spectator interactions. We 
observe that for negative values sin la as preferred by the global analysis of the unitarity 



triangle |5^, a measurement of the coefficient S'^tt could be used to determine 

sin 2a with a theoretical uncertainty of about ±0.1. Interestingly, for such values of sin 2a 
the "penguin pollution" effect enhances the value of the mixing-induced CP asymmetry, 
yielding values of S'tttt between —0.5 and —1. Such a large asymmetry should be relatively 
easy to observe experimentally. 

Although it illustrates nicely the effect of "penguin pollution" on the determination 
of sin 2a, Figure |Tl| is not the most appropriate way to display the constraint on the 
unitarity triangle implied by a measurement of 8^^.,^. In general, there is a four-fold 
discrete ambiguity in the determination of sin 2a, which we have reduced to a two-fold 
ambiguity by assuming that the triangle lies in the upper half-plane. Next, and more 
importantly, we have assumed that iKib/V^bl and /3 are known with precision, whereas a 
is undetermined. However, in the Standard Model |V^uf,/Kfe| and the angles a, /3, 7 are all 
functions of the Wolfenstein parameters p and f]. It is thus more appropriate to represent 
the constraint implied by a measurement of S'tttt as a band in the (p, f\) plane. To this 
end, we write 

_ PTin _2i^ ^ (1 - pf -fj^- 2ir/(l - p) ^ r^^ e^'^- 
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Figure 12: Allowed regions in the {p,f]) plane corresponding to constant 
values of the mixing-induced asymmetry 5"^,^, assuming the Standard Model 
(left), and using a fixed value sin20d = 0.48 (and cos20d > 0) for the B^-B^ 
mixing phase (right). The widths of the bands reflect the total theoretical 
uncertainty. The corresponding bands for positive values of S"^^ are obtained 
by a reflection about the p axis. The light circled area in the left-hand plot 
shows the allowed region obtained from the standard global fit of the unitarity 
triangle [Q. 



where r^^r parameterizes the large fraction in ( P^ without the prefactor 1/ Rh and is 
thus independent of p and fj. We now insert these relations into (pO|) and draw contours 
of constant S'^^r in the (p, fj) plane. The result is shown by the bands in the left-hand plot 
in Figure [1^. The widths of the bands reflect the total theoretical uncertainty (including 
power corrections). For clarity we show only bands for negative values of 5*7^77; those 
corresponding to positive S.„t, values can be obtained by a reflection about the p axis 
(i.e., f] —>■ — ff). Note that even a rough measurement of Sj^.,^ would translate into a rather 
narrow band in the (p, ?]) plane (this was also noted in |Q), which intersects the ring 
representing the |V^6/V^fe| constraint at almost right angle. This would therefore provide 
a very powerful constraint on the Wolfenstein parameters. We also stress that the sign 
of S'tttt, when combined with a measurement of I7I (using, e.g., the method described in 
Section |5.ij| ), can potentially determine whether the unitarity triangle lies in the upper or 
lower half-plane, and thus provide a nontrivial test of the CKM model of CP violation. 

The strategy outlined above remains useful even in the hypothetical case where Bd- 
Bd mixing is affected by new physics beyond the Standard Model (this scenario has 
recently be discussed in |5^, |60|, ^). Then the factor e~^*^ in the expression for 

Att^ in (^) must be replaced with the mixing phase e"^*'^'*, where (pd ^ P due to the 
presence of new physics. It is the value of sin 2(f)d that is measured in the time-dependent 
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CP asymmetry in the decays B^,B^ J/ipKs- From (pO|), we then obtain 

Aj^; = [ ± Vl - sin220, - ^ sin 20,] ^ ~ + "^"^^ , (92) 

which still implies a constraint in the (p, r^) plane for each pair of experimental values 
sm2(j)d. The ± sign refers to the sign of cos2(j)d. The result obtained for 
sin20d = 0.48 is shown in the right-hand plot in Figure |12[ In the figure we assume 
that cos20d > 0; the resulting bands for cos20rf < are, once again, obtained by a 
reflection about the p axis. We observe that the bands in the first quadrant of the (p, f]) 
plane intersect the rings from the iKit/V^fel constraint at almost the same places as in 
the Standard Model case, implying that for S^n < 1 (and cos20(i > 0) potential new 
physics effects in B^-Bd mixing have a minor impact on the results. On the other hand, 
the impact would be significant if S',r7r turned out to be positive. 



5.5 Predictions for CP-averaged branching ratios 

Whereas so far in this section we have focused on methods that require minimal input 
from the QCD factorization approach, we now discuss in detail our theoretical predictions 
for the B — > ttK, vrvr branching ratios. These predictions follow from the theory described 
in this work without relying on further phenomeno logical input. All branching fractions 
discussed in this section are averaged over CP-conjugate modes (even though for neutral 
K and B mesons this is not indicated by the notation). We use tb+ = 1.65 ps and 
tbo = 1.56 ps for the 5-meson lifetimes. 



Absolute predictions for branching fractions 



Two out of the seven decay modes, B"^ vr^vr^ and B^ -k'^K^ , are (almost) indepen- 
dent of the CKM phase 7, since the corresponding decay amplitude have to a very good 
approximation only a single weak phase. The predicted branching fractions for these 
modes are (setting 7 = 55° for concreteness) 



10^ Br(5^ 
10^Br(5^ ^ 7r^K° 
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0.28 



(93) 



where the first error is due to parameter variations as shown in Table |^, whereas the 
second one accounts for the uncertainty due to power corrections from weak annihilation 
and twist-3 hard spectator contributions. The dominant contributions to the uncertainty 
are shown in parentheses. Note that the decays B"^ tt^tt^ do not receive weak an- 
nihilation contributions. The B^ n^K^ branching ratio is to a good approximation 
proportional to (m^/llO MeV)~^'^^. The largest uncertainty for the n^n^ final state is 
a 50% normalization uncertainty from the current errors on \Vub\ and the B n form 
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factor. The systematics of theoretical errors is different for the t^K final states. Their 
absolute branching fractions are (approximately) proportional to the square of the pen- 
guin amplitude a^^ + r^a^ + TAh'i, which is sensitive to the penguin annihilation coefficient 
63 (see the left-hand plot in Figure The dominant theoretical errors therefore come 
from the annihilation parameter Xa and from the strange-quark mass (through r^). The 
theoretical errors detailed here are characteristic for the absolute branching fractions of 
all TTTT and ttK final states, respectively. 

The central values in (^) agree well with the current data summarized in Table 0. In 
particular, we find that the QCD factorization approach prefers large branching fractions 
for the ttK final states. This arises due to an enhancement of the QCD penguin amplitude 
relative to naive factorization. Weak annihilation contributions play only a minor role 
in this enhancement. Without annihilation the tt^K^ branching fraction would 

be reduced to 12.1 x 10^^. Hence the effect is not negligible, but there is no need for a 
largely enhanced annihilation contribution (this is in contrast to the findings of fl^])- In 
fact, the good agreement of our prediction with the data provides circumstantial evidence 
against the idea that annihilation effects could be much enhanced with respect to our 
estimates. For instance, increasing the parameter qa in (^) from 1 to 2 would increase 
the corresponding error on the branching ratio from to in which case it would 

require considerable fine-tuning of the strong-interaction phase of Xa to reproduce the 
experimental value of the branching ratio. 

It has also been suggested in the literature that one needs a large enhancement from 
penguin diagrams with a charm loop (so-called "charming penguins" ) to understand the 
overall nK branching fractions [Q. The leading perturbative contribution from the 
penguin diagrams in Figure ^, however, turns out to be very small. There exists also a 
power-suppressed contribution from these diagrams, when one of the quarks the gluon 
decays into becomes soft. In this case, the gluon is "semi-hard" and probes the charm- 
quark loop at a scale of order /i^ ~ (Aqcd "^fe)^''^- The penguin function G{s,x) in (^) 
tends to a constant for small x, which implies that the contribution from the semi-hard 
region to the function Gk{s) in ( p^D is suppressed by two (not one!) powers of the heavy- 
quark mass relative to the leading perturbative contribution. Hence, although the strong 
coupling constant in the semi-hard region is larger than 0^(7715), a large nonperturbative 
enhancement of the charm-penguin contribution appears implausible. (It is possible to 
obtain a first-order power correction by invoking a higher Fock component in the wave 
function of the emission meson, as discussed in [0. However, in this case the penguin 
loop continues to be a hard subprocess, and so contributes a factor of ^^(mft).) In a 
recent article a non-perturbative charming penguin contribution has been fitted to 
experimental data under the assumption that this effect is responsible for any deviation 
of the measured branching fractions from those expected within the Standard Model 
with p and f] determined by the standard unitarity-triangle fit. It it worth noting that 
such a fit is technically equivalent to fitting the annihilation contribution to the QCD 
penguin amplitude, indirectly related to Xa [eq. (0)] in our notation. If a modification 
of this amplitude were required (for which the present data does not provide strong 
motivation, as will become more evident below), we would attribute its physical origin 
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to weak annihilation rather than charm penguins, given the power-counting detailed 
above. 



In Figure |I3|, we show the dependence of the absolute branching fractions of the 
various B nK, nn decay modes (except for tt'^tt^) on the weak phase 7. In 

each plot, the solid line gives the central prediction of the QCD factorization approach 
at next-to- leading order in a^. For comparison, the short-dashed line shows the result at 
leading order, corresponding to naive factorization. The dark-shaded band is obtained by 
varying all input parameters as specified in Table ^, and by varying the renormalization 
scale between mb/2 and 2mb. It also includes the uncertainties due to the errors on the 
CKM parameters |Kb/14fe| = 0.085 ± 0.017 and \Vcb\ = 0.041 ± 0.003. The variation due 
to |V^fe/Kfe| alone is indicated by the long-dashed lines. The light-shaded band adds to 
this the uncertainties inherent to our modelling of power corrections due to twist-3 hard 
spectator and weak annihilation corrections, as discussed in Sections [4.2| and [4.3| . The 
annihilation contributions dominate the uncertainty in all cases (see Figure They 
imply a considerable uncertainty in the overall normalization of the irK modes. For the 
purposes of our discussion here the different sources of theoretical uncertainty are added 
in quadrature. Later, when our focus is on constraining CKM parameters, we will be 
more conservative and scan over the entire theory parameter space (corresponding to 
linear addition of errors). 

The next-to-leading order effects included in the QCD factorization approach signif- 
icantly enhance the branching fractions for the B — > irK modes with respect to their 
values obtained in the naive factorization model. No such enhancement occurs for the 
decays B^ tt+tt^ and B"^ tt'^tt^. As a result, the empirical finding that B — > ttK 
branching ratios are larger than the B ^ titi branching ratios is reproduced in our 
approach without any tuning of parameters. 

Predictions for ratios of branching fractions 

The uncertainty in the overall magnitude of the penguin amplitude due to weak annihila- 
tion effects, and the overall scale uncertainty in the predictions for the branching fractions 
due to hadronic form factors and CKM elements, are largely eliminated by taking ratios 
of branching fractions. The dependence of the six relevant ratios on the weak phase 7 is 
displayed in Figure in which the curves and bands have the same interpretation as 
in the previous figure. (The results presented here are consistent with the preliminary 
results published in [0, which did not include annihilation contributions. However, the 
treatment of theoretical uncertainties in differs from the present work.) Table |^ gives 



numerical predictions and theoretical uncertainties for some selected values of the angle 
7. It is evident that annihilation effects are less important for the ratios. The small error 
on the ratios involving nK final states shows the potential of these ratios to constrain 
7. Comparing our results with the preliminary experimental data collected in Tables |^ 
and H (and shown by the horizontal bands), we note a preference for large values of 7; 
however, the experimental errors are still too large to assign much significance to this 
observation. 
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Figure 13: Absolute CP-averaged branching fractions as functions of 7. Cen- 
tral values are shown by the solid line. The inner (dark) band corresponds to 
the variation of the theory input parameters (including iKib/KfeD? whereas the 
outer (light) band also includes the uncertainty from weak annihilation and 
twist-3 hard spectator contributions (see text for details). The long-dashed 
lines indicate the sensitivity to |Kife/V^6|. The short-dashed curve shows the 
result obtained using naive factorization. 
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Figure 14: Ratios of CP-averaged branching fractions as functions of 7. The 
meaning of the curves and bands is the same as in Figure [l^. The horizontal 
bands show the la and 2a error bands for the current experimental results. 
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Table 8: Predicted ratios of CP-averaged branching fractions for selected val- 
ues of 7. The last column shows the experimental values deduced from Table |^. 
Our averages ignore correlations between the individual measurements. 



Rat in 




n ."VT^l f^T*! TTl fiTl "I" 
J_/A.lJCi lliidilj 


2Br{7r"A'±) 

Br(7r±i^") 


0.94 ±0.07 1.16 ±0.07 1.44 ±0.16 1.70 ± 0.25 


1.41 ±0.29 


Br(7rTX±) 
2Br(7r"X") 


0.92 ±0.08 1.17 ±0.08 1.50 ±0.19 1.83 ± 0.34 


0.83 ±0.22 


Tgo Br(7r±i^0) 


0.74 ±0.07 0.91 ±0.04 1.12 ± 0.07 1.32 ±0.12 


1.06 ±0.18 


Br(7r+7r-) 
Br{7r=FA'±) 


0.96 ±0.60 0.67 ±0.38 0.43 ± 0.25 0.27 ±0.18 


0.26 ±0.06 


rg+ Br(7r+7r-) 
Tgo 2Br{7r±7rO) 


0.96 ±0.25 0.83 ±0.20 0.66 ±0.15 0.50 ±0.13 


0.42 ±0.14 



An important remark in this context concerns the ratio involving the neutral mode 
—>■ TT^K^, which seems to prefer a small value of 7. Our theoretical results indicate 
that to a good approximation (compare the two plots in the first row in Figure |l^) 

2Br(fi± ^ n^K^) _ Br(fi° ^ tt^K^) 

This relation is in slight disagreement (albeit by no more than 2a) with present data. 
Some authors have interpreted this fact as an indication of large rescattering phases 
(see, e.g., |6^). However, the validity of ( plf ) is a model- independent consequence of 
isospin symmetry, which is valid to linear order in the small tree-to-penguin ratios Q]. 
Therefore, we expect that with more precise data the values of the two ratios in ( P^ will 
come closer to each other. 

Finally, we wish to stress that the difference between the QCD factorization results 
and naive factorization are largest for the ratio of the tt+tt^ and tt^K^ final states. 
Whereas we can accommodate the low experimental value of this ratio for 7 > 50° (at 
the 2a level), this would not be possible for any value of 7 if one were to use the naive 
factorization model. 



It is evident from Figure that in many cases the error on |\4fe/V^f,| constitutes 
one of the largest uncertainties in the predictions for the ratios of branching fractions. 
This is true, in particular, for the interesting ratio of the tt+tt" and vr^i^^ final states. 
For instance, if |Kjb/V^{,| is 20% lower than its current central value, then our prediction 
for the 7t~^7t~ branching ratio is reduced by 40%, whereas the result for the 

B^ — i> n^K'^ branching fraction remains almost unaffected. This shows that for future 
analyses of rare B decays it will be important to reduce the theoretical uncertainty in 
the value of \Vub\- 

We can unfold the strong dependence of the ratios on |V^tif,/Kfe| by presenting the 
predictions for the branching fractions as functions of the Wolfenstein parameters p and 
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Figure 15: Two ratios of CP-averaged branching fractions as functions of the 
Wolfenstein parameters p and fj. The upper row corresponds to our central 
values; the middle (lower) row corresponds to the upper (lower) theoretical 
value including all uncertainties. We only show the upper half-plane (f/ > 
0). The current standard best-fit unitarity triangle and the annulus given by 
iKib/Kbl are overlaid to guide the eye. 
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f]. This is illustrated in Figure ^ for the two cases of the ratios 2Br(7r'^i^^)/Br(7r^A'°) 
and Br(7r+7r^)/Br(7r^i^^). The other two ratios involving irK final states exhibit a 
dependence similar to the first of these two examples. For not too large values of f/, 
any of the three ratios involving only ttK final states is a direct measure of p, since the 
dominant dependence on the CKM parameters arises from eKMCOS7 = tan^^c P- For 
positive p the theoretical uncertainty is relatively small, as can be deduced by compar- 
ing the three panels of Figure |1^ that refer to 2Bi(7r^ K"^) /Bt{7i'^ K^) . (Note that the 
contour lines in this plot do not coincide with those of constant inferred from Fig- 
ure since there the parameter £3/2 is assumed to be fixed and extracted from data, 
whereas here it is determined from theory and hence proportional to |V^6/V"cfe|.) The 
ratio Br(7r"'"7r~)/Br(7r^i^^) exhibits a rather different behaviour. The circular contours 
reflect the increase of the tt+tt" branching fraction with \ Vub\- The offset of the center to 
negative values of p results from the interference of the two contributions with different 
weak phases to the decay amplitudes, which for p > is constructive (destructive) for 
the TTTT (tiK) final states (and vice versa for p < 0). 

The current experimental values for the CP-averaged branching fractions are still 
afflicted by large errors, so a fit of [p, ff) to these ratios may appear premature. Nonethe- 
less, it is useful to carry out such an analysis at the present stage, not only to exhibit 
its potential once the errors decrease, but also to gauge the limitations eventually set by 
the theoretical uncertainties. Figure |l^ shows the results separately for each of the five 
ratios of CP-averaged branching fractions that have been measured so far. Consider first 
the dark-shaded band in each plot, which covers the region allowed by theory given the 
current central experimental value. The light-shaded bands in Figure |1^ represent the 
theoretically allowed regions obtained by using experimental values for the ratios shifted 
up or downwards by two standard deviations. Hence, the width of each band thus reflects 
the (largely irreducible) theory uncertainty given a certain value of the observable, while 
the current constraint on p and fj provided by each observable separately can roughly be 
taken to be the entire region bounded by the two hght-shaded bands. 



Global fit in the (p, 77) plane 

The information from all CP-averaged B nK, tttt branching fractions (and, in the 
future, from the corresponding CP asymmetries) can be combined into a single global 
fit, giving regions for the Wolfenstein parameters p and f] that are allowed by theory. 
We will now determine these regions, using two slightly different treatments of theo- 
retical uncertainties. Because our focus is to determine fundamental parameters of the 
Standard Model, we adopt a more conservative error analysis than in the previous para- 
graphs of this section and scan all theory input parameters over their respective ranges 
of uncertainty (rather than adding theoretical uncertainties in quadrature). The same 
conservative treatment was used in the analysis of the weak phase 7 in Section |5.3| . 

It is evident from our previous discussion that ratios of CP-averaged branching frac- 
tions provide the most powerful constraints on the Wolfenstein parameters, since they are 
afflicted with relatively small theoretical uncertainties. However, there are two reasons 
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Figure 16: Regions in the (p, fj) plane allowed by theory given the current 
central experimental values of the ratios of branching fractions (dark bands), 
and the experimental values ±2cr standard deviations (light bands). Only the 
upper half-plane {fj > 0) is shown. The dashed circles indicate the allowed 
region implied by the measurement of iV^ub/Kbl in semileptonic B decays. The 
light circled area shows the allowed region obtained from the standard global 



fit of the unitarity triangle |^ 
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Figure 17: 95% confidence region in the (p, f/) plane obtained from a fit to 
the CP-averaged B — >■ nK, nir branching fractions. The dark curve gives the 
envelope of the 95% confidence level contours, the lighter curve the envelope 
of the set of Xmm points. Here the solid lines refer to a cut on the xLin value, 
which corresponds to a 5% confidence level of the fit for a given model. The 
dashed lines refer to all models which have a Xmm P^^ degree of freedom of less 
than 1. See text for further explanations. 



for not performing a global fit to the ratios directly (see also the discussion in Appendix C 
of |5^): first, with the present large experimental errors there are significant differences 



in the results of the fit depending on whether the ratios or inverse ratios are used, lead- 
ing to an element of arbitrariness; secondly, even if two quantities have gaussian errors, 
their ratio does not. It is therefore preferable to perform the global fit for the individual 
branching ratios, despite the fact that the predictions for the branching fractions have 
larger (and correlated) theoretical uncertainties. 

We now present two versions of such a fit. In the first method, we find the allowed 
region from p and fj by determining first the 95% confidence level contour for a given 
theoretical input, and then scanning over all models in the theory parameter space that 
survive a 5% confidence cut. The theory parameter space is given by the set of all input 
parameters confined to the ranges specified in Table ^, variation of the renormalization 
scale II between mb/2 and 2mb, and scanning of the parameters Xh and Xa in the 
ranges described in (|7lD and (^). The resulting 95% confidence region is bounded by 
the dark solid lines in Figure |1^. The grey solid line inside this region envelopes the 
set of Xmin points of all models that pass the 5% confidence cut. The dashed dark and 
grey lines have the same interpretation except that only theory models with a Xmm P^^ 
degree of freedom of less than 1 are selected. The (p, f/) range so determined is compared 
with the standard CKM fit (also at 95% confidence level) which uses information 
from semileptonic B decays {\Vub\ and \ Vcb\), K-K mixing (e^^), and B B mixing {Attib^, 
AniB^, and sin 2/5). The constraint corresponding to iKi^/Kbl = 0.085±0.017 is indicated 
by the solid annulus. 

A fit to CP-averaged branching fractions cannot distinguish between f] and —f], and so 
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Figure 18: 95% (solid), 90% (dashed) and 68% (short-dashed) confidence 
level contours in the (p, ff) plane obtained from a global fit to the CP-averaged 
B ttK, tttt branching fractions, using the scanning method of [53]. The dark 
dot shows the overall best fit, whereas the light dot indicates the best fit for 
the default parameter set. See text for further explanations. 



for simplicity only the upper half of the (p, fj) plane is shown in the figure. At present our 
contours are symmetric about the p axis. In the future, data on direct CP asymmetries 
in the various B nK, tttt modes should be included in the fit, in which case it will 
become possible to obtain information about the sign of fj. 

Before interpreting the results of this analysis, we discuss a second method for per- 
forming the global fit, which follows the strategy outlined in [^]. It differs from the 
method explained above in two ways. First, the large overall scale uncertainty (of al- 
most 40%) due to the error on the value of the form factor F^^'"{0), which is common 
to the theoretical predictions for all branching fractions (see, e.g., (|93D), is eliminated 
by fitting the form factor to the data. In that way, the sensitivity to the Wolfenstein 
parameters in increased in the same way as it would be by using ratios; however, the 
statistical problems of using ratios (see above) are avoided. Specifically, we start from 
the function 

' Ei — S Tj 



X\S) = Y.[^^) , (95) 

where Ei are the experimental values for the branching ratios, cxj are their errors, and 
Ti are the corresponding theoretical predictions obtained with the fixed form factor 
F(f^'^(0) = 0.28. The scale factor S parameterizes the deviations of the product 
(F(f^'^(0) iKbl)^ from its central value. For a given set of theory parameters, we de- 
termine 

such that 'X^{S) is minimized, provided that the value of S is allowed by the theory 
ranges for the form factor (see Table 0) and \ Vcb\- For most choices of theory parameters 
this condition is satisfied. If it is not, we choose the value of S inside the allowed range 
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that is closest to the optimal value in (]96|). Next, we scan over all theory parameters 
and determine the global minimum of in the (p, ff) plane. We then plot contours of 
constant A^^ = ~ Xmin corresponding to the 68%, 90% and 95% confidence levels. 
The second difference with respect to our first scanning method is that these contours 
refer to a constant x^, whereas in the first method the is different for each of the 
theory models considered. The results of this analysis are shown in Figure 0. The best 
fit has an excellent = 2.1 (for 3 degrees of freedom), corresponding to a confidence 
level of 54%. The parameters of the particular theory model corresponding to this 
fit are all very reasonable. In particular, the quantity Xa parameterizing the weak 
annihilation terms is even smaller in magnitude than our default value. Hence, a good 
global fit to the data can be obtained without using extreme choices of input parameters, 
or invoking phenomenological recipes such as largely enhanced annihilation or charm- 
penguin contributions. The default parameter set in Table || also yields a good fit, 
which has = 4.57 and 21% confidence level. Note that at 90% confidence level the 
allowed range obtained by combining our results with the measurement of |V^f,/Kfe| in 
semileptonic B decays excludes = 0, thus establishing the existence of a CP-violating 
phase in 6 — > M transitions. 

The results of the two scanning methods shown in Figures and are similar 
(though the range obtained in the first method is somewhat more conservative). The 
allowed region for p and f\ obtained from our analysis of rare hadronic B decays is 
compatible with the standard global fit using information from semileptonic B decays, 
K-K mixing and B-B mixing. However, the best fits to the rare-decay branching 
fractions prefer a larger value of 7 (of about 90°). This has been noted by many other 
authors in the past, usually based on the naive factorization approximation and without 
a theoretical error estimate. Here we have put this analysis on a firmer theoretical 
footing. On the other hand, we see that small values of 7 are also allowed, provided 
the value of \Vub\ is lower than its standard central value. Combining our results with 
the standard fit would reduce the allowed region in the (p, ff) plane by about a factor 2, 
indicating that even with present experimental (and theoretical) errors the information 
obtained from rare hadronic decays provides important constraints on the Wolfenstein 
parameters. 

It will be interesting to follow how the comparison of the standard fit and the fit 
to rare hadronic B decays will develop as the data become more precise. If the two 
fits remain consistent with each other even as the allowed regions shrink in size, this 
would constitute a highly nontrivial test of the CKM model, in which for the first time 
the phase of the Vub matrix element (which enters through tree-penguin interference) 
is restricted by both direct and indirect measurements. If, on the other hand, the two 
regions would not overlap at a reasonable confidence level, this may suggest that the 
standard theory of weak interactions does not account completely for either the B-B 
mixing amplitude (since this is what determines the left-most side of the allowed region 
in the standard fit) or the loop-induced flavour-changing amplitudes in B ^ tcK, tctt 
decays. (Another option would be to abandon the theoretical framework advocated 
here, but we understandably leave it to others to pursue this avenue.) 
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5.6 Predictions for CP asymmetries 



The QCD factorization approach predicts the strong-interaction phases and hence the 
direct CP asymmetries in the heavy-quark hmit. In the following, we define the asym- 
metries as 

^^^^^^ - Br(i? ^ /) + Br(fi ^ /) ' ^^^^ 

where B = B^,B~ contains a h quark (rather than antiquark). The dependence of the 
CP asymmetries on the weak phase 7 is displayed in Figure |T9|. The asymmetries are 
typically predicted to be small, concurrent with the fact that strong-interaction phases 
are suppressed in the heavy-quark limit [l^. Even the largest asymmetries for the ttK 



final states are predicted to be of order 10% or less. An exception to this rule is the 
final state / = vr^vr^, for which the asymmetry is sensitive to QCD penguins and to 
the coefficient 02, both of which have large uncertainties and a potentially large relative 
phase. 

While our approach predicts the generic magnitude of the CP asymmetries, their 
precise values remain uncertain. This is not surprising, since in contrast to the branching 
ratios the asymmetries are more sensitive to - and in fact generated by - corrections to 
naive factorization. Hence, they are subject to larger relative uncertainties. In particular, 
the CP asymmetries are proportional to the sines of strong-interaction phases, which 
are of order asirrib) or Aqcd/'^6- Only the leading perturbative contributions to these 
phases are calculable. On the contrary, the CP-averaged branching fractions depend on 
the cosines of strong- interaction phases, which are equal to 1 in the heavy-quark limit. 

First measurements of CP asymmetries have been published by the CLEO |^ and 
Belle Collaborations. All results are compatible with no asymmetry, with typical 
la errors of about 20%. These errors are not yet small enough to draw meaningful 
conclusions, except that very large asymmetries appear already improbable. In the 
future, accurate measurements of CP asymmetries will test the generic prediction of the 
QCD factorization approach that strong-interaction phases are small. If this prediction 
is confirmed, further analysis of the CP asymmetries should help to determine the sign 
of 7, which cannot be probed using CP-averaged branching fractions. In addition, data 
on CP asymmetries may help to further constrain annihilation contributions and other 
theoretical input parameters in the QCD factorization approach. 



6 Comparison with other work 

A general overview of the various approaches to the evaluation of the matrix elements in 
exclusive hadronic B decays has been given in . In that reference we have commented 
on the strategies and methods that had been employed previously and explained how 
they are related to QCD factorization. In this section we discuss recent papers on the 
subject and clarify the differences with our work. We divide these papers into four 
categories: those using our QCD factorization formalism, those using a perturbative 
PQCD approach, those in which nonperturbative effects are estimated using light-cone 
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Figure 19: Direct CP asymmetries as functions of 7, assuming 7 > 0. For 
negative 7, the signs of the asymmetries are reversed. The meaning of the 
curves and bands is the same as in Figure |I^. The asymmetries vanish in 
naive factorization, so no corresponding hne is drawn. 



66 



sum rules or model calculations, and finally those containing phenomenological analyses. 
We consider each of these in turn. 

Before discussing other approaches to nonleptonic -B-decays it may be helpful to 
summarize again the conceptual basis of the QCD factorization approach. The key in- 
gredient is the systematic analysis of Feynman graphs in the heavy-quark limit, from 
which we deduce the factorization of infrared singularities into hadronic light-cone dis- 
tribution amplitudes and form factors. This enables us to establish the factorization 
formula (pSj). The separation of long- and short-distance contributions to the decay am- 
plitudes, necessary to establish factorization, holds only to leading power in Aqcd/'^6 
and is based on considerations analogous to those used to demonstrate factorization in 
other applications of QCD to hard processes (such as deep inelastic scattering, Drell-Yan 
production, and electromagnetic form factors of hadrons at large momentum transfer). 
The factorization formula leads to a model-independent treatment of exclusive hadronic 
B decays in the heavy-quark limit. A consistent counting scheme for powers of Aqcd/""^?) 
and a systematic identification of all the leading contributions are crucial for establishing 



this result. The framework proposed in M, O is general and provides a starting point 



for further theoretical developments, such as the improved understanding of the nonper- 
turbative input (e.g., the B ^ n form factor and the light-cone distribution amplitudes) 
and estimates of power corrections. 

6.1 Analyses within QCD factorization 

Several authors have used the framework of QCD factorization for applications to two- 
body hadronic B decays p|, |0> HI HI- We now compare the results reported 



in these papers with ours and comment on apparent differences and discrepancies. 

Muta et al. have generalized the results of our previous work 0] by including elec- 
troweak penguin contributions, and have applied the QCD factorization approach to the 
decays B nK, rnr. We would like to point out the following differences with respect 
to the present work: 

• The hard-scattering kernels are derived only for symmetric light-cone distribution 
amplitudes. 

• In evaluating the as corrections to the penguin coefficients 04 and ae, the existence 
of the two distinct types of penguin contractions shown in Figure § is not taken 
into account. As discussed in Section this leads to incorrect terms proportional 



to as{C4 + Cq) in these coefficients. 

An incomplete projector for the twist-3, two-particle distribution amplitude of the 
pion is employed. This gives an incorrect contribution proportional to in ag. 

We also observe a disagreement with the remaining terms in the correction of order 
as to ag, which concerns the function G'^.j^{sq) in [|I^. All three components, the Sq- 



dependent part, the constant, and the coefficient of the logarithm, differ from our 
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findings. We note that the //-dependence of G'j^^ is in conflict with the requirement 
of renormalization-group invariance of the product (fx) aQ^fi) . In addition, the 
contributions from the radiative vertex corrections to and ag are missing. 

A minor difference comes from the fact that we neglect electroweak penguin con- 



tributions to 04 and Og, while these are included in |]T8|. On the other hand, Muta 

"YcfT 
^77 



et al. omit the contributions proportional to Ci, C2 and Cj!^ to the electromagnetic 



penguin coefficients as and aiQ. 

The twist-3 contributions to the spectator hard-scattering amplitudes and annihila- 
tion effects are not discussed. As we have seen, these corrections are the potentially 
most important source of theoretical uncertainty. 



The formulae given in |]18| have also been applied to discuss B decays into vector- 



pseudoscalar final states |jT9[ and final states containing rj and rj' ||2^. The previous 
comments apply also to these papers, as they do to which relies on the same 
expressions except for dropping the formally power-suppressed terms proportional to ag 
and as- 

In a series of papers, Du et al. have discussed B decays into two light pseudoscalar 
mesons |2^. We focus on the two last papers in p2|, which contain the most complete and 



updated results. Similarly to the present work, electroweak penguin contributions and 
chirally-enhanced twist-3 components of the pion distribution amplitude are included 
in these papers, but no weak annihilation effects are considered. The hard-scattering 
kernels are given only for symmetric distribution amplitudes, and explicit results are 
presented for the case oi B ^ titt decays. We will therefore not distinguish between 
TTTT and ttK final states in the present discussion. The formulae for the coefficients 
ttj can be directly compared with our results. A minor difference concerns the tiny 
electroweak penguin contributions to the coefficients ai,...,aQ, which we decided to 



neglect in our approximation scheme, but which are retained in |Q . Next we recall from 
Section that in the approximation of including twist-3 contributions only when they 
are chirally enhanced, the equations of motion require the use of asymptotic distribution 
amplitudes $p, Du et al. have considered a more general situation, which would 
affect the cancellation of endpoint singularities and the renormalization-scale dependence 
of some results. Use of the appropriate asymptotic distribution amplitudes eliminates 
these spurious ambiguities. Finally, we note two minor discrepancies: 



Our expression for the twist-3, hard spectator interaction in (|57|) contains a factor 

(98) 



2/ip ^Mix) ^p{y) 



rrib X y 

in the integrand, whereas the corresponding result given in eq. (37) of of the third 



paper in has 

Mb X Qf ' ^ ^ 
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Using the asymptotic forms of the twist-3 distribution amphtudes $p and $0-, we 
find that the two results differ by a factor of y. 



In our resuh for the twist-3 penguin contribution P3 in ([5^), the coefficient of Ci, 
for example, is 

5 '■'(?) 4 -'^(^''' « 



The equivalent expression given in eq. (23) of the third paper in ||22| reads 



+ ^ A.) 1^ (^) + ^ + ^ - ^'(^^) - ^'^(^p) • (101) 

Using asymptotic twist-3 distribution amplitudes, one may check that G'{sp) + 
G"{sp) = G{sp) and A^^ = 1/2. The two results then agree up to a constant 1/6. 
This difference can be traced back to an inconsistent use of twist-3 projectors in 
four space-time dimensions within a rf- dimensional loop calculation, before the sub- 
traction of ultraviolet poles is performed (see the last reference in |^2[). Analogous 
differences occur in the other terms in and ag. 

Apart from these discrepancies, the expressions agree with our results for symmetric 
distribution amplitudes. 

Cheng and Yang have applied the QCD factorization approach to a study of the 
decays B —>■ (pK [^. Annihilation topologies are discussed and argued to be important 



numerically. A few minor discrepancies with our results occur in the expressions for 
the coefficients a3,...,io. The penguin contraction of the operators (^4,6 is treated incor- 
rectly, as discussed in Section |37^ . Photonic penguin contractions of the operators (5i,2 
contributing to ag 10 are omitted. Also, QCD corrections to ag^g are not considered, pre- 
sumably because the corresponding amplitudes are formally power suppressed. Similarly 
to the results of ||22|, the twist-3 kernel for the hard spectator interaction contains an 
additional factor of y. Finally, we briefly comment on the analysis of weak annihilation 



contributions. In ||23| the final state consists of a vector meson and a kaon rather than two 
pseudoscalar mesons. The decay amplitudes are estimated in the approximation of using 
leading-twist distribution amplitudes for the meson, but including chirally-enhanced 
twist-3 contributions for the kaon. Taking this systematic difference into account, we 



can still compare the results of with our expressions. The two agree in the case of 
{V — A)®(y — A) operators, but the results for {S — P)®{S + P) operators are different. 
The latter are the ones that depend on the twist-3 contributions, but we do not have 
sufficient information to trace back the origin of the discrepancy. 



6.2 The PQCD approach 

Recently, Keum, Li and Sanda have presented an extensive study of B irK decays 
within a perturbative hard-scattering (or "PQCD") approach ||T^. While the underlying 
goal of a separation of soft and hard physics in the i?-decay matrix elements is similar 
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in spirit to QCD factorization, there are fundamental differences in tlie implementation 
of this idea. 

In the PQCD approach, the B ^ n and B ^ K form factors are assumed to 
be perturbatively calculable. This assumption is justified by invoking Sudakov effects 
to regulate the infrared-sensitive contribution from the endpoint region in the integral 
over the light-cone momentum of the outgoing spectator quark. In other words, the 
contribution from the region where the spectator quark is soft is supposed to be strongly 
suppressed, and therefore the exchange of a hard gluon is always required. Instead of 
being of leading order (in the QCD coupling constant), as in the generic case of a i? — vr 
form factor dominated by soft physics, the form factor is now counted as being of order 
Us in perturbation theory. Thus, the hierarchy of the various contributions to the decay 
amplitudes in PQCD is very different from that in our approach: 

• In the PQCD approach, the B irK matrix element is a quantity of order a.^- To 
this order, therefore, all the Wilson coefficients can then be taken in the leading 
logarithmic approximation. 

• The "nonfactorizable" hard gluon-exchange contributions to the kernels T/ in (p5|) 
enter at order and are therefore omitted. Note, however, that it is essential to 
consider such effects in order to establish factorization. This has never been done 
explicitly in the PQCD framework. 

• The form-factor terms containing the kernels T/ and the hard-scattering terms 
containing the kernels T/^ in the factorization formula are treated as being of the 
same order in the PQCD approach. The former are called "factorizable" and 
the latter "nonfactorizable". Because both terms vanish simultaneously in the 
limit as —>■ 0, naive factorization is not recovered in any limit of QCD. This is in 
contrast to our interpretation of the QCD factorization formula, for which naive 
factorization is recovered in the heavy-quark limit. 

After this brief synopsis of the key features of the PQCD scheme, we now examine 
critically the main assumptions of this approach in the context of B —> tiK decays. 



Importance of Sudakov effects and calculability of the S — > tt form factor 

In the PQCD approach the transverse momenta of the quarks are kept explicitly when 
evaluating contributions that are potentially infrared sensitive. For instance, the gluon- 
exchange kernel for the B ^ tt form factor is written as 

H ^ ^ ^ + . . . , (102) 

{x^ml + /c|_L)[xiX3mg + {ki± - hx^] 

where the ellipses represent less singular terms. Here Xi (xs) is the light-cone momentum 
fraction of the spectator quark in the B meson (pion). In the conventional collinear 
expansion adopted in our work, the transverse momenta k'^j_ ~ Aqqj^ would be treated 
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as higher-twist effects and dropped (more precisely, the amphtude would be expanded in 
powers of transverse momenta). It then follows that H ~ 1/3:3, and so the convolution 
with the leading-twist pion distribution amplitude $7r(x3) ~ ^3(1 — 0:3) would lead to an 
integral / dx^/x^, which is infrared divergent. This well-known result can be interpreted 
as a sign of the dominance of soft endpoint contributions in the B ^ tt form factor 



63. 

In [|17| (see also |^) this divergence is regulated by keeping the transverse momenta. 
Then a Fourier transform from transverse momentum space into impact parameter space 
(6 space) is performed, and a Sudakov factor is included for each of the meson distribution 
amplitudes. This factor strongly suppresses the region of large b, and the B ^ ir form 
factor is therefore assumed to be perturbatively calculable. 

It should be noted that retaining the transverse momentum dependence in H, which 
is only part of the complete higher-twist contribution, is a model-dependent procedure. 
It is precisely in the critical endpoint region X3 ^ that the leading two-particle Fock- 
state description that underlies this analysis breaks down. A further puzzle to us is the 
assumption made in that the spectator quark in the B meson has a large momentum 
component in the "— " direction, fcf = Xim^/v^, when the pion momentum is in the 
"+" direction. This is used to derive the Sudakov factor for the B meson in analogy to 
the case of a fast light meson. Since the spectator quark in the B meson is intrinsically 
soft, we see no justification for such an assumption. 

In conclusion, we believe that the perturbative evaluation of the B ^ it form factor. 



which is one of the central ingredients of the PQCD analysis in is not justified. 
The relevance of Sudakov form factors in hadronic B decays should be investigated 
in a more systematic way. However, it seems unlikely to us that Sudakov logarithms 
are sufficiently important at the scale mb to eliminate the soft contributions to heavy- 
to-light form factors and thus render them calculable in a model-independent way. A 
complete description of the B —>■ tt transitions at large recoil in the heavy-quark limit 
and the derivation of a (hypothetical) factorization theorem for these processes has, to 
our knowledge, never been presented. For recent progress in this direction and further 
critical discussions of the problem see [p3, BO]. 



Dynamical enhancement of penguin contributions ("fat penguins") 

Keum et al. claim a dynamical enhancement of (factorizable) contributions to the matrix 
elements of penguin operators (such as (^4 and Qq) in the effective weak Hamiltonian. 
The first reason given in support of this assertion is a dynamical enhancement of the 
B ^ 71 form factor of the scalar density ub contained in (the Fierz-transformed form of) 
Qq, relative to the form factor of the vector current uj^b relevant for other operators. 
A dynamically different structure for the form factor of the scalar density in PQCD is 
claimed, and a strong sensitivity of the effect to the model employed for the i?-meson 
light-cone distribution amplitude is noted. However, an enhancement of this nature is 
not possible, because the form factor of the scalar density is related to that of the vector 
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current by the equations of motion, leading to 



2 



(n+lublB^) = ^ . (103) 

nib 



The erroneous conclusion in is the consequence of an incorrect treatment of the 
twist-3 contribution to the pion distribution amplitude (see below). 

The second reason given in favor of a penguin enhancement is the choice of a low 
scale fi in the (leading-order) penguin coefficients a^lfi) and aQ{fi). This is motivated 
by arguing that /i should be a typical scale intrinsic to the dynamics of the B ^ n 
form factor in the PQCD evaluation. However, the scale n is unphysical and must be 
canceled by vertex corrections to the operator matrix elements, and in the case of a6(/x) 
also by the running quark masses. Such corrections have been properly included in the 



present work, but they were neglected in |Tj]. The scale n in 04^6 may depend on the 



scale of non-factorizable hard spectator interactions, but it is clearly independent of the 
internal dynamics of the form factor. The arguments made in favour of a low scale n 
and correspondingly increased penguin coefficients are therefore not well justified. 



Relevance of annihilation topologies 

In their analysis, Keum et al. find that penguin annihilation contributions are numeri- 
cally very important and give a dominant contribution to some of the B ttK decay 
amplitudes. We have investigated these effects in the present work and found them 
to give a significant correction (compatible with being a power correction of canoni- 
cal size), but not a dominant contribution. Nevertheless, it must be emphasized that 
weak-annihilation diagrams contribute at subleading power in the heavy-quark expansion 
and, in general, are not calculable in a QCD-based factorization approach. Therefore, 
the question about the relevance of these effects warrants further investigation. 



Generation of large, calculable strong-interaction phases 

It is claimed that large strong-interaction phases are generated by hard gluon exchange 
with the spectator quark in the B meson, as well as by weak annihilation diagrams. Let 
us consider the hard spectator interactions in detail. In that case the effect is calculable 
within the QCD factorization approach and found to be real to leading order. The 
source of the imaginary part found in is as follows. The quark propagator entering 
the hard-scattering diagram is written as 



Working to leading power we would drop the transverse components and find a real 
contribution to the kernel. Keum et al., on the other hand, keep the transverse momenta 
and hence generate an imaginary part proportional to 5(fc5_ — X2X^m1), where fc_L = 
/c2± + k^±- Again, this treatment is model dependent. The Fourier transform of the 
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kernel into impact parameter space results in a Bessel function with the imaginary part 
proportional to ruhb jQ{y/x2X^mbb) . In the evaluation of the matrix element this function 
is convoluted with the Sudakov factor. Since the Bessel function is oscillating rapidly, 
with the amplitude of m^fo Jq growing like v^mfefe, the result of this convolution is very 
sensitive to the details of the 6-space cut-off. This implies that the estimate of the strong- 
interaction phase is both model dependent and numerically sensitive to effects that are 
poorly under control. Similar comments hold for the estimate of the strong-interaction 
phases from the annihilation contributions, which are generated in an analogous way. 

Finally, we remark that the simple 75 structure of the twist-3 projection for the 
pion employed by Keum et al. is incomplete. The proper treatment is discussed in 
Section p75[ The wrong twist-3 projection is, in particular, inconsistent with gauge 
invariance. Moreover, the correct asymptotic behaviour of the twist-3 pion distribution 
amplitude ^p{x) is proportional to a constant, whereas the functional form ~ x(l — x) 
is assumed in |T^. This problem affects all decay amplitudes, including the B ^ it 
form factor, the corresponding spurious penguin enhancement, and the nonfactorizable 
spectator interactions and annihilation contributions. 



Other works using the PQCD approach 



Other analyses in the spirit of the PQCD approach were presented in |]70|, |72|, [T^ , 
to which similar comments apply. In ||72|, the presence of a "recoil-phase" effect was 



advocated, which was claimed to originate within the PQCD framework. This phase 
should affect, e.g., the B ^ it form factor, which was assumed to be dominated by 
hard gluon exchange. It was argued that, when the gluon is exchanged between the 
spectator quark and the b quark, the 6-quark propagator could go on the mass shell 
because niB > mt + ?Tispect- This in turn would lead to a complex phase in the form 
factor. In our opinion this conclusion is unwarranted, since bound-state effects have to 
be factorized before a perturbative treatment can be justified and thus cannot influence 
the hard-scattering process. Another fundamental objection to the "recoil phase" is that 
it would contradict the fact that the B ^ n form factor is a real quantity. 



6.3 Estimates of nonperturbative effects 

Recently, Khodjamirian has suggested to study hadronic matrix elements for two-body 
B decays in the framework of light-cone QCD sum rules |^. As an example, he discusses 
the matrix elements of the current-current operators Qi and Q2 for B — > vr+vr" decays. 
The idea is to generalize the light-cone QCD sum-rule analysis of the B ^ 71 form factor 
directly to the case of hadronic two-body modes. The starting point is the correlation 
function 

Faip, q,k) = -J d'x e^^P-")^ J A e'^P-'^y (0| T [jj^\y)Qmj'^''\m^-i<l)) , (105) 

where jj^^ = U'jaJbd and j^^^ = mhbi'j^d are interpolating currents for the emission 
pion and the B meson, respectively. The explicit pion state vr~(g) represents the recoil 
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pion that absorbs the spectator quark. According to the QCD sum rule philosophy, the 
correlator is evaluated in two ways: by a direct calculation in QCD, and by inserting 
complete sets of hadronic states between the operator and the interpolating currents, 
extracting the desired ground-state contribution with the help of quark-hadron duality. 
The two sides of the equation are expressed in the form of dispersion relations, and the 
usual Borel transformation is applied. To leading order in a^, the factorized result for the 
matrix element is recovered. A particular type of power correction is then estimated as 
a further illustration. This contribution comes from higher twist, three-particle (quark- 
antiquark-gluon) Fock components of the recoil pion, where the (nonfactorizable) gluon 
couples to the quark lines of the emission current. In the heavy-quark limit the resulting 
expression is demonstrated to scale as AQCD/''^b relative to the leading contribution. 
It is estimated to be of relative size XeItub with A^; ~ 0.05-0.15 GeV. Note that this 
correction has no rescattering phase. Its numerical effect is small, but comparable to 
the perturbative corrections at leading power. Additional nonfactorizable contributions 
exist but have not yet been investigated. Examples are the gluon-exchange effects that 
correspond to the vertex corrections and the spectator scattering diagrams in the QCD 
factorization formula. In the sum-rule approach these effects are in principle contained 
in similar diagrams, together with certain hadronic corrections of subleading power. A 
potential difficulty will be to disentangle power corrections to the asymptotic result given 
by the QCD factorization formula from uncertainties intrinsic to the sum rule method, 
such as the assumption of quark-hadron duality and the approximation of the emission 
pion by an interpolating current. It is not fully clear how this can be achieved in a 
controlled and systematic fashion. 

Some recent papers |7^, ^] have tried to model soft final-state interactions via 
the rescattering of certain hadronic channels such as B ^ DD vrvr. We consider 
such an approach to be problematical for a variety of reasons. There are many more 
intermediate channels beyond those taken into account. Systematic cancellations among 
these channels, which are predicted to occur in the heavy-quark limit, are missed when 
only a few intermediate states are retained (see the discussion in Sections 3.4 and 7.2 of 
p!5|). Moreover, the hadronic dynamics of multi-body decays is very complicated and in 
general not under theoretical control. The main problem of such models is the lack of a 
systematic approximation scheme based on parametric expansions. In our opinion, the 
use of a purely hadronic language, suitable for kaon decays, is not very helpful in the 
case of B decays, where the number of channels and the energy release are large. 

In [Q, the coefficients a, obtained from the QCD factorization approach have been 



used in conjunction with a hadronic description of final-state interactions. The general 
caveats concerning hadronic rescattering models apply also here. Moreover, such an 
approach faces a manifest double-counting problem, since rescattering effects in the 
heavy-quark limit are already contained in the QCD coefficients a^. 

A nonperturbative treatment of "charming penguins", matrix elements with charm 
loops of the type shown in Figure ^ was proposed in |7^. In this calculation operators 
with flavour structure (c6)(sc) were split into nonlocal products of currents (c6) and (sc). 
These were connected by intermediate virtual D-meson states to yield transitions such 
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as -B — s> -D — s> Kt: under the separate action of the currents. However, the operators 
(c6)(sc) are local operators in the effective theory at the 6-mass scale. They would 
become nonlocal only if probed at the far higher scale of Mw-, but certainly not at a 
scale of order m^^, as implied in We therefore see no basis for the method adopted 
in this paper. 

6.4 Phenomenological analyses 

There are several recent analyses in the literature that approach hadronic B decays in a 
more phenomenological way. Some of the most extensive studies of this type have been 
presented by Ali et al. |]TB[, who elaborate on the concept of "generalized factorization" 
|79| . In these analyses, part of the order-a^ vertex and penguin contributions calculated 



in the present work are included. As a consequence, the decay amplitudes receive strong- 
interaction phases due to the Bander-Silverman-Soni mechanism ||3^, and nonzero CP 
asymmetries are obtained. However, no systematic attempt to include all such effects (or 
to prove factorization) is made. Also, these authors rely on ad hoc model assumptions 
such as an "effective number of colours" N^. ^ 3, which is introduced to parameterize 
nonfactorizable corrections. 

A recent study of S — tiK decays presented in emphasizes general, model- 



independent parameterizations of the decay amplitudes such as our parameterization in 
(p!8|), avoiding as much as possible the use of dynamical input. For instance, strong- 
interaction phases are a priori allowed to take any value in the construction of bounds 
on CKM parameters. Information on these phases may then be inferred indirectly. 
Such an approach is complementary to the more ambitious goal of exploiting theoretical 
insight into the hadronic matrix elements. In fact, the smallness of the rescattering 
effects parameterized by e*"^" in (|18D , which is a prediction of the QCD factorization 
approach, represents a useful input to such a phenomenological analysis. 

Many other phenomenological studies of two-body hadronic B decays have been per- 
formed in the literature. Recent examples are |81|, |8^, p^, p3, pR p6| , where simplifying 



assumptions are made and no systematic estimates of theoretical uncertainties are un- 
dertaken. 



7 Conclusions 

In this paper we have presented a detailed study of -B — > t^K and B ^ Tin decays based 
on the QCD factorization formula. This approach allows us to perform a systematic and 
model- independent calculation of two-body hadronic B decays in the heavy-quark limit. 
We have evaluated the hard-scattering kernels entering the matrix elements for B nK, 
7171 decays at next-to- leading order in ag, thus obtaining predictions for the decay am- 
phtudes including the leading, "nonfactorizable" corrections. We have included the con- 
tributions from electroweak penguins, taking into account the corresponding order-a^ 
corrections to the dominant terms (enhanced by mf/M^r or 1/ sin^O^r) in a consistent 
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approximation scheme. All hard-scattering kernels have been derived for general, asym- 
metric distribution amplitudes, as appropriate for K mesons. 

In addition to computing the model-independent leading-twist results, we have iden- 
tified and estimated those power corrections which are expected to be the largest. We 
have analyzed the complete set of contributions from light-cone wave functions of twist 3 
with a chiral enhancement factor / {mu+md) or mj^ / {mg+nis) , as well as the contribu- 
tions from weak annihilation topologies, including both twist-2 and twist-3 components 
in the light-meson wave functions. We distinguish three classes of power corrections: 
the contributions from operators with (pscudo-) scalar currents (~ oe.s), twist-3 effects 
in the hard spectator interactions, and weak annihilation corrections. While the first 
two effects are proportional to the chiral enhancement factor, the annihilation terms are 
quadratic polynomials in this factor. Since these are power-suppressed effects, we do 
not expect the factorization formula to work in these cases. In fact, naively evaluating 
these terms in a "hard-scattering" framework one encounters logarithmic endpoint sin- 
gularities. (Incidentally, the terms of order related to ag.s a-re free of such infrared 
singularities, but divergent terms are present for the other two contributions.) There- 
fore, while our results at leading twist are model-independent predictions of QCD in the 
heavy-quark limit, some model dependence is currently unavoidable in the description of 
power corrections. We choose to regulate these divergences by introducing a simple cut- 
off, leading to the complex phenomenological parameters Xh and Xa- The motivation 
for this model-dependent approach to power corrections is twofold: we obtain systematic 
order-of-magnitude estimates, while automatically keeping track of relative suppressions 
or enhancements from Wilson coefficients and CKM factors. In addition, we can in- 
vestigate the sensitivity of the observables on the model parameters. It turns out that 
annihilation graphs may in principle give sizable corrections. On the other hand, twist-3 
effects in the hard spectator contributions appear generally to be less prominent. 

Using this framework, we have performed a detailed comparison with the available 
experimental data on B ^ ttK and B ^ tttt branching ratios. The data are beginning 
to be sufficiently precise to allow for detailed phenomenological analyses. An important 
result is that in the QCD factorization approach the B irK branching fractions can 
be larger than the B ^ nn ones without any tuning of the theoretical input parameters, 
and without invoking large phenomenological power corrections. An acceptable fit to the 
branching fractions is obtained even if we impose that 7 < 90° as implied by the standard 
constraints on the unitarity triangle. This is in contrast to the naive factorization model, 
in which the observed branching fractions can either not be reproduced at all or require 
a large CKM angle 7, in confiict with other indirect determinations. Encouraged by 
this success, we have derived constraints in the (p, f]) plane from a global fit to the 
B — > ttK, tttt branching ratios, which already provide useful additional information on 
the unitarity triangle. Clearly, the experimental situation will continue to improve and 
allow us to further test our framework and exploit it in phenomenological studies of 
hadronic B decays. An important role in these analyses will be played by the tree-to- 
penguin ratio £3/2, which can be reliably determined from B^ — > ir^ir^ and B^ — > tt^K^ 
decays and can be unambiguously confronted with the theoretical predictions. This will 
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be a valuable check that the power corrections are not surprisingly large. The absence of 
penguin and annihilation contributions in the B"^ — > tt^tt^ mode is a crucial feature in 
this analysis. A more precise measurement of this decay will thus be particularly useful. 

We would like to conclude by emphasizing again the strategy underlying our treat- 
ment of weak annihilation contributions and non-leading, but chirally-enhanced specta- 
tor interactions. While the present data do not require these effects to be larger than 
expected, the data also do not definitively exclude this possibility. Our error analysis is 
based on allowing the parameters qa and qh to be smaller than 1 (reflecting the range of 
our expectations), but the error estimates sometimes depend sensitively on the precise 
choice of this upper limit, in particular in the case of qa- With more precise data one 
may contemplate fitting qa to data (in practice this implies fitting the QCD penguin 
amphtude) in order to decide upon the plausibihty of the adopted range of values. The 
discovery of large annihilation contributions would by itself constitute valuable insight 
into the strong-interaction dynamics of nonleptonic decays, but would evidently limit the 
utility of the QCD factorization approach. The present data make this scenario appear 
unlikely, but we look forward to further information to validate our error estimation. In 
this case, important and reliable information on flavour physics will become available 
from B — > nK and S — > tttt decays, as well as from a large class of other two-body 
hadronic decays into light pseudoscalar and vector mesons. 
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